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SUMMARY 


The objective of this research is to develop models and methods to analyze 
the interaction between flight and propulsion systems. In this report, 
perturbation equations which describe flight dynamics and engine operation 
about a given operating point are combined to form an integrated 
aircraft/propulsion system model. The equations used to describe aircraft 
dynamics include dependence of aerodynamic coefficients upon atmospheric 
variables and, as a result, altitude is used as a state variable. The engine 
equations are derived from a low order engine model. In order to describe the 
interface between flight and engine variables, an off-design engine 
performance model is used to develop perturbation equations which describe the 
effect of flight condition and inlet performance upon propulsion system 
variables. These perturbation equations are used to identify interaction 
parameters in the integrated system model. 

Linear quadratic regulator methods and numerical linear algebra techniques 
have provided a flexible means to analyze interaction effects upon the 
stability and control of the integrated system model. The inclusion of 
interaction effects in the model changes the overall system behavior. To 
analyze interaction effects on control, consideration is first given to 
control requirements for the airframe and engine models separately. For the 
separate airframe model, feedback control provides substantial improvement in 
the short period damping. For the integrated airframe/propulsion system 
model, feedback control compensates for coupling present in the model and 
provides good overall system stability. Analysis of suboptimal control 
strategies indicates that performance of the closed loop integrated system can 
be maintained with a feedback matrix in which the number of non-zero gains is 
small, relative to the total number of components in the feedback matrix. A 
method based on eigenvalue sensitivity analysis has proved to be an effective 
means for determining which gains in the integrated system feedback matrix can 
be set to zero while, at the same time, maintaining system performance. 



INTRODUCTION 


In present and future aircraft designs consideration will be given to 
problems associated with thrust management, fuel efficiency, improved 
maneuverability and pilot workload . Concepts such as variable geometry 
inlets and engines along with thrust vectoring and reversing provide a degree 
of interaction between propulsive and aerodynamic forces that requires a more 
complete integration of airframe and propulsion control systems. Improved 
performance is also predicted through further integration of flight controls 
with guidance, fire control, weapon delivery, and structual control systems . 
This trend tov/ard integration of subsystem controls has motivated the 
development of models and methods to analyze the behavior of dynamically 
coupled subsystems and to design control laws for the improvement of subsystem 
cooperation and the enhancement of overall system performance. 

Flight propulsion interface . Airframe/propulsion interactions are a major 
concern in aircraft design and many complex problems are associated v/ith the 
description of this interface. Major problems are associated with the 
description of external airflow effects on the airframe, inlets, engines, and 
nozzles. Problems also occur in the development of aerodynamic accounting 
systems and in the verification of performance criteria. Many aspects of 
these problems have been previously considered, and references [1, 2, 3] 
indicate the scope of past activities. Two NASA programs which have dealt 
with interaction problems are the Cooperative Control Program (YF12) and the 
Integrated Propulsion Control Program (IPCS). The interaction between flight 
and propulsion system is broadly described by the effect of flight condition 
variables upon propulsion system mass flow, pressures and temperatures and by 
the effect of propulsion system forces and moments upon the aircraft. 
Examples of severe engine/inlet/airframe subsystem interactions have been 
observed in NASA flight research programs. An indication of the types of 
interactions observed is given in table 1. 

TABLE 1. Airframe/Propulsion interactions. 


1 r~ 

] Aircraft | 

1 1 

Interaction 

1 n 

[ Result 

Ref 


[ F-104 1 

airframe/inlet 

1 divergent lateral oscillations 

4 


1 F-lll 1 

engine/inlet 

1 distortion factor exceeds limits 

4 


! YF-12 1 

airframe/ inlet/ engine 

I unstable dutch roll and phugoid 

4 


1 F-15 1 

1 L_ 

airframe/ inlet 

1 improved static stability 

j _ ,j 

5 



Both open and closed loop interactions have been observed between airframe, 
inlet, and engine. These interactions are a consequence of airflow variations 
between engine and inlet, airframe forces and moments induced by the 
propulsion system, and variations of inlet/engine operation with flight 
condition. Subsystem interactions have had a direct effect on aircraft 
stability, control, and performance. Further aspects of interaction effects 
are discussed by Schweikhard and Berry [4] and by Hunt, Surber and Grant [6]. 
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Control configured aircraft with variable geometry engines and inlets will 
utilize interaction effects by design. In this application, consideration 
must be given to the physical interaction betv;een flight and propulsion state 
variables and the cross-coupling effect of subsystem controls. 

Models and control . The development of models of subsystem interactions for a 
given control requirement is complicated by the complexity of the system under 
consideration. In a simulation of the YF-12 aircraft [7] , six degree of 
freedom aircraft motions were represented along with a three axis stability 
augmentation system. The effect of mixed compression, variable geometry 
inlets on aircraft forces and moments was simulated along with effects of 
aircraft motions on the behavior of inlets in the started mode. Further 
effects of changes in flight condition upon engine operating conditions were 
also considered in the simulation. Other simulations and models of the 
flight/propulsion interface are discussed by Tinling and Cole [8] and by Cole, 
Sellers and Tinling [9] . 

A low order model of a variable geometry inlet and a turbofan engine were 
interfaced in the work of Michael and Farrar [10] . Optimal control methods 
were used to design a closed loop controller for the integrated system. The 
control of an integrated airframe/propulsion system by state regulation was 
previously considered by W. R. Seitz [11]. 

For interacting systems, problems arise in deciding which interactions to 
include in a model and in the development of appropriate models of the 
interface between subsystems. The necessity of considering an integrated 
model which incorporates a full set of propulsion system control parameters 
along with the airframe control parameters has been discussed by Sevich and 
Beattie [12] . Sevich and Beattie indicate that a manageable design approach 
is to first optimize propulsion system control based upon overall system 
requirements. Then add to the aircraft control parameters those engine 
control parameters which effect the given control requirements. Control of 
the integrated system model should then be considered. The decomposition of a 
control problem for an interacting system into problems for lower dimensional 
subsystems is a common approach. The subsequent solution of control problems 
for each subsystem and their combination into a solution for the overall 
system has been called the decomposition principle. Further examples and 
background on the decomposition principle appear in the work of D. D. Siljak 

[13] . 

The objective of this research is to develop models and methods to analyze 
the interaction between flight and propulsion systems. In this report, 
perturbation equations are obtained which describe flight dynamics and engine 
operation about a given operating point. A model of the standard atmosphere 

[14] is used to describe changes in ambient temperature, pressure, air density 
and speed of sound with altitude. The equations used to describe aircraft 
dynamics include dependence of aerodynamic coefficients upon atmospheric 
variables and, as a result, altitude is used as a state variable. The engine 
equations are derived from a low order engine model, which was used in control 
design studies by DeHoff , Hall, Adams and Gupta [15] . In order to describe 
the interface between flight and engine variables, an off-design engine 
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performance model is used to develop perturbation equations which describe the 
effect of flight condition and inlet performance upon propulsion system 
variables. These equations are used to develop an operating point model of an 
airframe/propulsion system at a given flight condition. Although a single 
integrated airframe/propulsion system model is analyzed in this report, the 
methods used to obtain this model apply at other flight conditions. 

Linear quadratic regulator methods and multivariable system analysis 
techniques are used to analyze subsystem interaction effects on stability and 
control. For the given operating point, feedback controllers are designed for 
the separate flight and engine models. These results are used to aid in the 
design of a feedback controller for the integrated airframe/propulsion model. 
Feedback controllers for the integrated system designed by linear quadratic 
regulator methods involve a large number of nonzero gains and consequently, 
result in a complex control system with many active controls. In this report, 
various stategies to obtain less complex control laws are compared. The 
analysis indicates that the inclusion of subsystem interactions can change the 
overall stability and control of the system. Analysis of suboptimal control 
strategies indicates that control system performance can be maintained using 
feedback matrices with a small number of nonzero gains. 

Acknowledgement . The principal investigator. Dr. Robert E. Fennell, spent a 
sabbatical leave from Clemson University at NASA Langley Research Center in 
1979-80. During this period the principal investigator worked in close 
cooperation with Mr. F. J. Lallman of NASA Langley Research Center in the 
development of the problem and methods of analysis discussed in this report. 
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SYMBOLS 


Air. Ctr. Aircraft with feedback control. 
Aircraft 

perturbation variables 

V velocity (m/sec) 

a angle of attack (rad) 

q pitch rate (rad/sec) 

0 pitch attitude (rad) 

h altitude (km) 

y glide path angle (rad) 

6e horizontal stabilator (rad) 

M Mach number 

nominal variables 

Vq velocity (m/sec) 

Mg Mach number 

hg altitude (km) 

cl, c2, c3 weighting coefficients 

D /Dt differentiation with respect to time 

E( ) expected value 

Eng. Ctr. engine with feedback control 
Engine 

perturbation variables 


N1 

fan speed (rpm) 



N2 

compressor speed (rpm) 



P5 

augmentor pressure (kPa) 



Wf 

main burner fuel flow (N/sec) 


P2 

compressor discharge pressure 

(kPa) 

Wfc 

command fuel flow (N/sec) 



A 

nozzle area (sq.m) 



CIW 

inlet guide vane (deg) 



RCW 

rear compressor variable vane 

(deg) 

BLC 

compressor bleed (%) 



Th 

net thrust per engine (N) 



Wa 

fan airflow (kg/sec) 



T4 

turbine inlet temperature 

(K) 


SMAF 

fan stall margin 



SMHF 

compressor stall margin 



DP/P 

relative fan exit pressure 

change 

Pr 

inlet pressure recovery ratio 
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SYMBOLS 


nominal variables 

TH thrust (N) 

WA fan airflow (kg/sec) 

WF/WA fuel to air ratio 

ST specfic thrust (thrust per unit mass flow) 

r^ j inlet pressure recovery ratio 
/ 

GSM(g) gain significance matrix 

Im. imaginary part of complex number 

Int. open loop integrated system 

Int. Ctr. integrated system with feedback control 

J, JA, JE, JI performance indices 

K feedback gain matrix 

Re. real part of complex number 

Sen(X,e) relative sensitivity of X with respect to e 
Stl, St2, St3, St4 suboptimal control strategies 
T modal matrix 

th time to half amplitude 

tp time of period of oscillation 

ua, ue, u aircraft, engine and integrated system control vectors 

xa, xe, X aircraft, engine and integrated system state vectors 

xO initial vector 

Xi fundamental mode 

ya, ye, y aircraft, engine and integrated system response vectors 
GREEK SYMBOLS 

6. pressure at engine location i relative to sea level 

^ standard (101.4 kPa) 

E system parameter 

0. temperature at engine location i relative to sea level 

^ standard (288.2 K) 
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SYMBOLS 


A real Jordan canonical form 
X eigenvalue 

5 damping coefficient 

wn natural frequency (rad/sec) 

SUBSCRIPTS 

i engine location 

s static 

SUPERSCRIPTS 

1 first row of a matrix 

' matrix transpose 


INTEGRATED SYSTEMS 


Control system design problems arise in the analysis of large scale systems 
composed of many interconnected subsystems. Use of an integrated system model 
allows the investigation of interaction effects in the control design process. 
In this manner, physical interactions between subsystem state variables and 
control variable cross-coupling effects can be considered from the outset in 
the design process. 

Design methods frequently focus upon operating point conditions, which 
provide a simplified system model. Linearization procedures and model 
reduction techniques provide further simplifications. The interaction between 
a given subsystem and other subsystems at a given operating point may be 
described by the linear system of equations 

D(XI)/Dt = AI»XI + BI*UI + E_,_ CIJ«YJ + DI*VI 

J rl 

YI = FI*XI + GI«UI + HIJ«YJ + EI*ZI 

J r I 

where XI, UI, VI, ZI denote the state, control, and noise variables for the 
I-th subsystem and the YJ denote interaction variables between the subsystems. 
Here AI , BI , CIJ, DI , FI, GI , HIJ , and El denote appropriate system matrices. 
The matrices CIJ and HIJ are referred to as interconnection matrices in this 
report and describe the interaction between subsystems at a given operating 
point . 

In order to introduce the coordination of subsystem controls into the 
design process an adequate model of subsystem interactions is necessary. 
Frequently these interactions are difficult to describe and the determination 
of adequate models of interaction effects requires detailed analysis of 
subsystem properties. The problem of determining parameters in an interaction 
matrix from given system properties and data is commonly referred to as a 
parameter identification problem. A method to determine subsystem interaction 
matrices for an airframe/propulsion system model will be presented in this 
report. Details of this method are presented in appendix A. In this example, 
interaction matrices are determined in a mannner so that the steady state 
responses of the system to changes in interconnection variables agree with 
previously developed steady state design point models. 

In this report, linear quadratic regulator methods and multivariable system 
analysis techniques are used to analyze the effect of subsystem interactions 
upon overall system stability and control. Combination of subsystem models 
into an integrated system model results in a system in the standard form 

Dx/Dt = A*x + B*u + D*v 
y = F»x + G»u + E«z. 

Control systems of this form have received extensive study. The stability 
and response of the unperturbed system is determined by the eigen-structure 
(eigenvalues and eigenvectors) of the system matrix A. The representation of 
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solutions of this system is described in most system analysis texts [16], 
[17]. The following approach is used in this report. The real Jordan 
canonical form of A is denoted by A. The matrix T denotes the modal matrix 
such that 


A*T = T«A. (1) 

Solutions of the initial value problem 

Dx/Dt = A*x x(0) = xO 


may be written as 

x(t) = e^^»xO 

= T*e^^«T“^»xO 

= cl*Xl(t) + c2*X2(t) + ... +cn»Xn(t) 

where 


T(cl, c2, ... ,cn)‘ = xO 


and 


T*e^^ = (Xl(t), X2(t), ... ,Xn(t)). 

The functions Xi(t), i=l,..., n, are the fundamental modes of the system and 
the structure of the modal matrix T determines the interdependence between 
system components and the fundamental modes, and consequently the 
interdependence between system components. Coupling between subsystems caused 
by the inclusion of subsystem interactions will be illustrated in the 
airframe/propulsion system model which follows. 

Frequently subsystem interactions have been treated as small perturbations 
and neglected in control design studies. The sensitivity of system 
eigenvalues and eigenvectors to small perturbations has been extensively 
studied, an introduction to this analysis is given in the text by G. W. 
Stewart [18] . Small perturbations can drastically change the system eigen- 
structure. Equation (1) may be used to analyze the sensitivity of the eigen- 
structure to changes in system parameters. If e denotes a system parameter 
then 


(DA/De)T + A(DT/De) ^ (DT/De)A + T(DA/De) 

whenever these derivatives exist. These equations may be solved to determine 
eigenvalue and eigenvector senitivities , details are presented in references 
[19], [20]. 

For the case of distinct eigenvalues the canonical form K, the modal 
matrix T, and the fundamental modes Xi, i = l,...,n, have a particularly 
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simple structure. In this case changes in the system eigen-structure caused 
by the inclusion of interaction matrices may be readily calculated and the 
determination of DA/De is straightforward. For completeness these details are 
summarized in appendix B and an identity is presented which avoids complex 
arithmetic in the calculation of the sensitivity of complex eigenvalues. 

Sensitivity analysis techniques may be used to study interaction effects 
upon open and closed loop stability and control. In this report the effect of 
a system parameter upon stability will be described by the relative changes in 
system eigenvalues due to changes in the system parameter. Thus, if X and e 
denote a system eigenvalue and parameter respectively then the sensitivity of 
X with respect to e is measured by 

Sen(X,£) = (DX/Ds)*le/X| . 

To a first approximation, one obtains 

AX = (DX/De)Ae 


and consequently 


AX/IXl = Sen(X,s) [As/ I e I ] . 

Thus Sen(X,e) is a measure of the relative change in X due to a change in e. 

Linear quadratic regulator methods [16] , [17] will be used to determine 

stable regulators for interacting systems. Systems will be written in the 
standard form 


Dx/Dt = A»x + B»u 
y = F»x + G*u 

with a performance index of the form 


J = /r \Y'*W-»y + u'«U «u dt. 

[0,*o)‘‘ 1 1 

The control which minimizes this performance index is 

u = -u”^(B' *5 + R' )x 

where Q satisfies the matrix Riccati equation 


A'«Q + Q*A + W - (Q*B+R)U“^(B‘ »Q+R' ) = 0 

or equivalently 

(A-B*U"^»R' ) 'Q + Q(A-B*U"^*R' ) + F ' (W^-W^»G«U~^«G ' *W^)F - Q»B*u"^»B*»Q = 0 
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where 


W = F'«W^*F 
R = F'»W^*G 
U = + G' -W^^G. 

In this case if x(0) = xO then J = xO’*Q«xO. In general if u = K«x is a 
stabilizing feedback control then J = xO'»Q«xO where Q satisfies the Liapounov 
equation 

(A+B»K)'Q + Q(A+B*K) + (F+G»K) ' (F+G»K) + K'*U^»K = 0 

In either case if xO is considered as a random initial condition with 
covariance matrix cov(xO,xO) = I then the expected value of J, denoted E(J), 
may be used as a scalar performance index [17, p. 371] and 

E(J) = Trace(Q). 

A numerical package of Fortran coded subroutines, ORACLS [20], was used to 
perform the numerical linear algebra required in this report and to solve 
linear quadratic regulator problems. The solution of the regulator problem 
uses principally one of two routines from the ORACLS package, RICTNWT and 
ASYMREG. RICTNWT uses the Newton-Kleinman algorithm to solve the continuous 
steady-state Riccati equation. Some numerical problems arose when RICTNWT was 
used to solve the regulator problem for the airframe/propulsion example of 
this report. The algorithm did not converge. As a result the routine ASYMREG 
was used to solve the associated Riccati equation and with this procedure 
satisfactory results were obtained. 
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INTEGRATED AIRFRAME/PROPULSION MODEL 

A description of an integrated airframe/propulsion system operating point 
model is presented in this section. All equations represent variations about 
the given operating point and all variables should be interpreted as 
perturbations about the given operating point. Longitudinal flight conditions 
are considered. 

Airframe . The linearized aircraft longitudinal equations of motion are of the 
form 


Dxa/Dt = AA*xa + BA«ua + CAE«ye 


ya = FA*xa 


The state, control, and response variables and typical equilibrium values for 
this model are as follows: 


xa. V, velocity (265.6 m/sec) 

xa2 a, angle of attack (0.0761 rad) 

xa^ q, pitch rate (0. rad/sec) 

xa. 0, pitch attitude (0.0761 rad) 

xa- h, altitude (13.72 km) 

D 

ua^ 6e, horizontal stabilator (-0.0346 rad) 

ya M, Mach number (0.9) 

ya^ h, altitude (13.72 km) 

ye^ Th, net thrust per engine(12 833 N) . 


The interaction matrix CAE describes the effect of thrust upon the 
longitudinal state variables. Data in this report is derived from a twin- 
engine, advanced fighter aircraft model in which it is assumed that thrust 
acts parallel to the aircraft centerline and the engines are located a short 
distance below the center of gravity. In the derivation of these equations a 
model of the standard atmosphere [14] was used to describe changes in ambient 
temperature, pressure, air density and speed of sound with altitude. The 
dependence of aerodynamic coefficients upon these variables resulted in the 
inclusion of altitude as a state variable. It will be noted later in this 
report that the resulting model is unstable. 


System matrices for the given operating point are listed in table 11 of 
appendix C. This data was supplied by NASA. 

Engine . Multivariable control design methods have been used extensively in 
the design of control laws for turbofan engines. The report [15] contains a 
discussion of current control design methods. In large scale problems, model 
reduction techniques may be used to provide simplified system models. The 
model studied in the present report is derived from a low order engine model 
which was obtained through model reduction techniques [15] and has the form 
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Dxe/Dt = AE*xe + BE»ue + CEA»ya + CEI*yi 
ye = FE*xe + GE*ue + HEA*ya + HEI»yi. 

The state, control, and response variables and typical equilibrium values for 
this model are as follows: 


xe- N1 , fan speed (9 785 rpm) 

xe^ N2, compressor speed (12 401 rpm) 

xe^ P5, au^entor pressure (74.89 kPa) 

xe^ Wf, main burner fuel flow (3,31 N/sec) 

xe^ P2, compressor discharge pressure (673.8 kPa) 

5 

ue. Wfc, command fuel flow (3.31 N/sec) 

^ 2 
ue^ A, nozzle area (0.259 m ) 

ue CIW, inlet guide vane (0.997 deg) 

ue^ RCVV, rear compressor variable vane (3.99 deg) 

ue_ BLC, compressor bleed (0.997%) 

ye. Th, net thrust (12 833 N) 

ye^ Wa, fan airflow (27.58 kg/sec) 

ye^ T4, turbine inlet temperature (1 590 K) 

ye^ SMAF, fan stall margin(0 . 14164) 

ye^ SMHC, compressor stall margin(0 . 1445) 

ye^ DP/P, relative fan exit pressure change (0.9973). 

Corresponding system matrices are listed in table 12 of appendix C. 

Inlet. In this example the aircraft inlet is characterized by a single 
parameter, Pr, inlet pressure recovery ratio. The presure recovery ratio is 
defined as the ratio of total pressure at the engine face to the total 
freestream pressure. Optimal engine performance requires high pressure 
recovery. The main requirement of an aircraft inlet is to provide mass-flow 
of sufficient uniformity to the engine to preclude stalls and to operate at 
high pressure recovery. In general, the subsonic inlet must supply air at the 
engine face at a specified Mach number and without separation of the flow. 
The supersonic inlet must decelerate the flow. Variable geometry inlets are 
currently designed to match variations in airflow demanded at the engine face 
with variations in flight condition and engine operating condition. Changes 
in inlet geometry affect not only engine performance but also aerodynamic 
forces and moments. Variable geometry inlet effects are not included in the 
model presented in this report. 

Integrated model. A simplified off-design performance model of a dry turbofan 
cycle de\^loped by F. J. Lallman [21] is used to analyze the effect of flight 
condition and inlet pressure recovery upon engine state and response 
variables. This off-design performance model expresses the engine specific 
thrust, fuel to air ratio, along with total temperatures and pressures 
throughout the engine, as a function of Mach number, ambient temperature, 
inlet pressure recovery ratio, fan pressure ratio, compressor pressure ratio. 
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and burner temperature. Specification of fan pressure ratio, compressor 
pressure ratio, and burner temperature yields a design point model from which 
engine performance and state variables are determined as functions of Mach 
number, altitude, and pressure recovery ratio. Details of this analysis are 
contained in appendix A.- 

For a turbofan engine with a fan pressure ratio of 2.9, a compressor 
pressure ratio of 7.93 and a burner temperature of 1559 K at a flight 
condition with a Mach number of 0.9 and altitude of 13.72 km, the design point 
model yields the following perturbation equations: 


P5 = 

42.020 

M - 

15.938 

h 

+ 

101.016 

Pr 

Wf = 

2.717 

M - 

0.5905 

h 

+ 

3.749 

Pr 

P2 = 

635.55 

M - 

90.631 

h 

+ 

574.812 

Pr 

Th = 

8 267 

M - 

2 121 

h 

+ 

17 826 

Pr 

Wa = 

25.678 

M “ 

4.35 

h 

+ 

27.586 

Pr 


The interaction matrices CEA, CEI, HEA, and HEX in the engine model are 
determined so that the steady state response of the variables augmentor 
pressure, P5, compressor pressure, P2, thrust per engine, Th, and fan airflow, 
Wa, to step changes in Mach number, M, altitude, h, and pressure recovery 
ratio, Pr, agrees with the above relationships. Details of this procedure are 
presented in appendix A and values for the interaction matrices for the given 
operating point are contained in table 12 of appendix C. It should be noted 
that the dependence of fuel flow, wf, upon flight condition has been removed 
in the determination of interaction matrices and consequently changes in the 
overall system behavior should be expected. 

Combination of the airframe and engine models leads to the following 
operating point model for the integrated system: 

Dx/Dt = AAAE»x + BABE*u 
y = FAFE*x + GAGE*U 

where 



AA + CAE -HEA* FA 

CAE*FE 


BA 

CAE*GE 

AAAE = 

CEA* FA 

AE 

BABE = 

0 

BE 









_ 


and 



FA 

0 


0 

0 

0 

FAFE = 

HEA 

FE 

GAGE = 

0 

GE 

HE I 




- 




_ 


For the given operating point values of the matrices appear in table 13 of 
appendix C. In this model the term CAE*HEA*FA is due to the dependence of 
thrust upon Mach number and altitude, the term CAE»FE results from the 
dependence of thrust on engine states, and the term CEA*FA results from the 
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dependence of engine temperatures and pressures on Mach number and altitude. 
The term CAE*GE represents changes in thrust due to engine controls. 

Interaction effects . The inclusion of subsystem interactions in the 
airframe/propulsion system model has changed the overall system behavior. 
This change in the system behavior is due to changes in the eigen-structure of 
the open loop system. Separate and integrated system eigenvalues are listed 
in table 2, below, whereas the separate and integrated system modal matrices 
are listed in table 14 of appendix C. Eigenvectors corresponding to the 
engine in the separate and integrated system models are essentially unchanged 
whereas eigenvectors corresponding to the airframe have changed significantly. 
The magnitude of the difference between corresponding eigenvectors of separate 
and integrated systems along with the angle between these vectors is listed in 
table 2. It should be noted that corresponding normalized eigenvectors are 
being compared, i.e. eigenvectors of unit length. 

TABLE 2. Interaction effects on stability. 


System eigenvalues. 


Airframe 


j Separate 
1 Integrated 
1 

-2.76E-3 

1.91E-3 

3.74E-4 ± j3.25E-2 
-3.65E-4 ± j3.65E-2 

1 -6.78E-1 ± 
1 -6.78E-1 ± 

j2.20E+0 

j2.20E+0 

1 

1 

1 

1 

1 

1 Separate 
1 Integrated 
1 

1 

-5.81E-1 

-5.63E-1 

Engine 

-1.88E+0 1 -6.59E+0 
-1.88E+0 ! -6.59E+0 

1 -l.OOE+1 i 
i -l.OOE+1 j 

-1.72E+2 

-1.72E+2 

1 

1 

1 

1 

1 

_| 

1 

1 System eigenvectors. 
1 

1 Distance 1.37 

1 Angle (deg) 86.4 

1 

Airframe 
1.37 I 1.37 

86.2 1 86.2 

1 1.28 1 

1 79.6 1 

1.28 

79.6 

n 

1 

1 

1 

1 

1 

1 

1 

1 Distance 

1 Angle (deg) 
1 

0.00 

0.00 

Engine 

0.00 1 0.00 

0.00 1 0.00 

1 0.00 1 

1 0.00 1 

0.00 

0.00 

1 

1 

1 

1 

J 

The significance of 
by a comparision of 

these changes in the system 
the transient response of 

eigen-structure is evidenced 
the separate and integrated 


models. The separate models do not include coupling between airframe and 
engine state variables and, consequently, changes in the states of one system 
do not effect the states of the other system. Also, the response of the open 
loop, integrated system to small offsets in engine state variables indicates 
weak coupling in the system. In this case, the open loop engine states 
respond almost exactly as in the separate model and small oscillations, 
relative to the given operating point values, are induced in the airframe 
states. This behavior is expected due to the structure of the system 
eigenvectors. Typical responses of the separate and integrated system models 
to an offset in fan speed, Nl, are depicted in figure 1. The response of the 
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model to small changes in airframe states illustrates the effect of coupling 
present in the integrated system model. Within a small time period, offsets 
in the values of velocity, v, and altitude, h, can lead to significant changes 
in engine state and output variables. Also the airframe response, in this 
case, differs slightly from that of separate model because thrust variation 
with altitude has been included in the model. Offsets in the values of angle 
of attack, a, pitch rate, q, and pitch attitude, 0, have little effect upon 
engine variables. The response of the separate and integrated models to 
offsets in values of angle of attack and altitude are depicted in figures 2 
and 3. Here a one per cent offset in the value of altitude leads to 
equivalent changes in the values of augmentor pressure, P5, compressor 
pressure, P2, thrust, Th, and inlet airflow, Wa. It should, also, be noted 
that within the given time period some variables have not attained their 
maximum displacement. 

In general, sensitivity analysis may be used to study the effect of 
subsystem interconnections upon system stability. If e denotes a parameter in 
the system matrix A(e), then the sensitivity of an eigenvalue \i to the 
parameter e can be measured as 

Sen(Xi,e) = DXI/De • Ie/XII. 


The parameter e may be a parameter in an interconnection matrix or some other 
interconnection parameter. For example, consider the system matrix 


A(s) = 


AA + E CAE»HEA«FA e»CAE»FE 


e*CEA«FA 


AE 


which is obtained by replacing CAE, CEA, and HEA by e»CAE, e«CEA, and e«HEA 
in the integrated system model. Here e =0 corresponds to the separated 
models while e = 1 corresponds to the integrated model. Calculations for the 
example model indicate that Sen(Xl), Sen(X2), and Sen(X3) are large in 
comparision to Sen(Xi), i=4, — ,10. Consequently, the first three modes are 
most sensitive to this interconnection parameter. These are the modes 
associated with the phugoid motion and the mode introduced by consideration of 
altitude as a state variable. Values of DXi/DE and Sen(Xi,£) for e = 1 are 
listed in table 3. 
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TABLE 3. Sensitivity to subsystem interactions. 



Xi 

"I 

DXi/De 1 

Sen(Xi,l) 

n 

1 

_l 

1.912E-3 



7.411E-3 ] 

3.875E+0 

n 

-3.654E-4 

+ 

j3.647E-2 1 

-5.175E-4 + j7.554E-3 j 

“1.416E+0 + j2.071E-l 


-3.654E-4 

- 

j3.647E-2 1 

"5.175E-3 - j7.554E-3 j 

-1.416E-0 - j2.071E-l 


-5.628E-1 



-2.137E-3 I 

“3.797E-3 


-1.883E+0 



3.771E-4 1 

2.002E-4 


-6.781E-1 

+ 

j2.200E+0 1 

9.565E-5 - j2.252E“4 j 

1.411E-4 -- jl.024E-4 


-6.781E-1 

- 

j2.200E+0 1 

9.565E-5 + j2.252E-4 j 

1.411E-4 + jl.024E-4 


-6.587E+0 



-3.896E-3 1 

-5.915E-4 


-l.OOOE+1 



5.543E-19 1 

5.543E-20 


-1.722E+2 


1 

-5.524E-3 I 

“3.208E-5 

J 


17 


SYSTEM CONTROL 


Linear quadratic regulator methods are used to design a feedback control 
law for the integrated airframe/propulsion system and subsystem interaction 
effects upon the resulting feedback control law are analyzed. Consideration 
is first given to the control requirements for each subsystem at compatible 
operating conditions. Feedback control laws are determined for each subsystem 
and these results are used to aid in the design of a feedback control of the 
integrated system. To allow full consideration of all interactions full state 
feedback is considered for the integrated system. Control of the integrated 
system involves several active controls which are determined by a feedback 
control law with many non-zero gains. Suboptimal control strategies are 
compared to determine if performance can be maintained with simplified 
feedback control laws. In particular, a method which uses sensitivity 
analysis is found to provide an effective means for determining simplified 
feedback control laws. 

Flight control . The objective of the flight control problem is to improve 
longitudinal performance through coordinated utilization of interacting 
aerodynamic and propulsive forces and moments. For the operating point of 
this example the short period damping is marginal. The performance index used 
is of the form 


JA = /p vya‘*WA*ya + ua'*UA*ua dt 

Lo,<»; 

where the response vector ya = FFA*xa + GGA»ua has components 

ya^ = V, velocity 

ya^ = y, glide path angle 

ya^ = q, pitch rate 

and the control vector has components 

ua^ = 6e, horizontal stabilator 
ua^ = Th, thrust per engine. 

Weighting matrices WA and UA were chosen with a small penalty on the use of 
thrust, in order to allow thrust control in the longitudinal model. The 
particular weighting matrices and other system matrices are listed in table 11 
of appendix C. 

Solution of the optimization problem provides the feedback control relation 
ua = -K*xa with 


K = 


8.11E-5 3.15E-1 -3.89E-1 -2.48E-1 -3.72E-4 

8.25E+0 -5.02E+3 5.40E+3 4.46E+3 3.14E+2 


and xa = (v, a, q, 0, h) . The large gains in this feedback matrix are due to 
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the choice of units for thrust. If thrust were measured in kilo-Newtons then 
each entry in the second row of K would be scaled by 0.001. 

The improvement in performance obtained through the use of this feedback 
control may be illustrated by comparing the stability, performance, and 
transient response of the system 

Dxa/Dt = (AA - E*BA*K)xa 

where e is a parameter which ranges between zero and one. The value e = 0 
corresponds to the open loop system while e = 1 corresponds to the closed loop 
system with feedback relation ua = -K*xa. The transition in eigenvalues and 
short period performance may be observed by varying e between zero and one. 
The eigenvalues of AA -s»BA*K with e = 0.0, 0.2, 0.4, 0.8, and 1.0 are listed 
in table 4. 

TABLE 4. Longitudinal eigenvalues AA - e»BA»K. 


E 1 

1 


1 

0.0 

1 -2.76E-3 

0.2 1 

-2.90E-3 

0.4 

-3.03E-3 

0.6 

[ -3.13E-3 

0.8 

1 -3.21E-3 

1.0 

1 -3.24E-3 


1 


Eigenvalues 


3.74E-4 ± j3.25E-2 
-1.88E-2 ± j2.65E“2 
-1.95E-2 , -5.50E-2 
-l.llE-2 , -9.95E-2 
-8.34E-3 , -1.38E-1 
-6.97E-3 , -1.77E-1 


T 

I -6.78E-1 ± j2.20E+0 
I -l.OOE-1 ± j2.14E+0 
I -1.22E+0 ± j2.02E+0 
I -l,49E+0 ± jl.85E+0 
I -1.76E+0 ± jl.61E+0 
I -2.03E+0 ± jl.26E+0 

j 


Calculation of the short period natural frequency, damping coefficient, 
period, and time to half amplitude illustrates the effect of the feedback 
control on longitudinal performance. Changes in short period performance are 
listed in table 5. The substantial improvement in short period damping should 
be noted. 

TABLE 5. Short period variables. 


wn - natural frequency (rad/sec) 
tp - period (sec) 

— 1 1 "1 


^ - damping coefficient 

th - time to half amplitude (sec) 


«n 


0.0 

0.2 

0.4 

0.6 

0.8 

1.0 


2.303 

2.336 

2.360 

2.377 

2.385 

2.385 


0.2945 

0.4054 

0.5154 

0.6254 

0.7364 

0.8496 


tp 


2.729 

2.690 

2.662 

2.644 

2.634 

2.635 


th 


1.022 
0.732 
0.570 
0 . 466 
0.395 
0.342 
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The improved damping and stability obtained from the feedback control is 
further illustrated by a comparison of the transient response of the closed 
loop system with that of the open loop system. The transient response of the 
open and closed loop models to an offset in each state variable is depicted in 
figures 4 through 8. Several properties of the closed loop system are 
apparent in these figures. The improved short period damping due to the 
weight on pitch rate, q, in the performance index is evidenced by the absence 
of oscillations in the closed loop system response. In each case, it appears 
that altitude, h, is being traded off for improvements in velocity, v, and 
glide path angle, y. Also, it appears that a small amount of thrust control, 
Th, is being used. Although it is not apparent from the figures, it should be 
noted that the closed loop system is stable and all variables will return to 
equilibrium. 

To determine the significance of thrust control in this example, a second 
regulator problem was posed with only the horizontal stabilator, 6e, as a 
control. Use of a performance index with WA as previously defined and UA = 4 
resulted in a feedback control law fie = -Kl»xa with 

K1 = [ 7.69E-5 3.17E-1 -3.89E-1 -2.49E-1 -5.42E-4 ] 

The gains in this feedback control are very similar to those of the first row 
of the previously determined feedback matrix, K. The eigenvalues for the 
closed loop system AA -BA»K1 are -2.24E-3, -6.97E-3, -1.77E-1 and -2.03±jl.26. 
These eigenvalues should be compared with those of the closed loop system AA - 
BA«K (see table 4). The only difference in eigenvalues is in the stability of 
the mode associated with the introduction of altitude as a state variable. 
This mode is most closely associated with the phugoid motion. It should also 
be noted that the improvement in short period damping previously obtained is 
also obtained by use of the control fie = -Kl«xa. Comparisons of the transient 
response of the two closed loop systems indicates that the two systems respond 
in a similar manner to offsets in state variables. The main difference 
between the two systems is that velocity responds faster in the system with 
thrust control, i.e. AA - BA*K, than in the system with the horizontal 
stabilator as the only control, AA - BA»K1. Typical responses of the two 
closed loop systems are compared in figure 9. 

Engine control . Detailed analysis of the regulator design for the engine is 
contained in reference [16]. At this operating point the engine is in a 
region of fan stability sensitivity to augmentor ignition. The performance 
index adapted from reference [16] is 

JE = / r xXe'«WE*xe + ue'*UE*ue dt 

Lo,oo) 

where xe and ue represent the state and control variables of the engine model 
(see p. 13). Table 14 of appendix C contains values for the matrices WE and 
UE. The resulting gain matrix is 
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K = 


9.02E-4 1.40E-3 6.58E-4 3.48E-1 

3.72E-5 2.03E-6 -8.58E-4 -1.36E-3 

-8.99E-3 3.33E-3 1.03E-1 -6.31E-1 

6.78E-4 -7.98E-3 2.21E-2 3.96E-1 

-4.02E-5 -5.86E-5 -1.32E-3 -3.17E-2 


-7.93E-4 

-1.78E-5 

-8.38E-3 

2.30E-2 

-6.58E-4 


Integrated system control . In order to pose a regulator problem for the 
integrated system a weighted performance index of the form 


J = cl«JA + c2»JE + c3*JI 


is used where JA, JE are as before and 

J1 - J r xYI ' -YI dt 

[O, oo) 

denote performance indices for the airframe, engine and inlet respectively and 
cl, c2 and c3 denote the weights. In the model example, YI denotes the single 
parameter Pr, inlet pressure recovery ratio. This performance index may be 
written in the standard form 

J = ,.y'*W*y +u'«U»u dt 

[o , oo; 

where y = FF*x + GG»u with y = (v, Z, q, Th, N1 , N2, P5 , Wf, P2) ' , x = (v, a, 
q, 0, h, Nl, N2, P5, Wf, P2) ' , u = (6e, Wf c , A, CIVV, RCVV, BLC, Pr) ' and 



1 

> 

o 

_J 


" 0 

0 

0 

FF = 

(HEA*FA)^ (FE)^ 

GG = 

0 

(GE)^ 

(HEI)^ 


0 I 


0 

0 

0 




and 



cl*WA 0 

0 “ 


~ cl*UA^^ 

0 

0 

w = 

0 CI.UA 22 

0 

u = 

0 

c2*UE 

0 


0 0 

c2*WE 


0 

0 

c3*I 


_ 




— 


— 


The superscript 1 in the above equations is used to denote the first rov; of 
the indicated matrices. In this model thrust is no longer a control variable. 
For each choice of weights a feedback control 

u = -K(cl, c2, c3)»x 

is determined. The feedback obtained by normalizing JA and JE about their 
previously obtained optimal values and by placing a large weight on JI (cl = 
.318, c2 = .338, c3 = 10E6) is 
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K = 


4.587D-02 1.7220+00 -6.M85D-01 -2.7300+00 -5. 4800+00 

1.364D-02 -9.8220-02 -6. 1780-03 5. 6800-02 5.1180-02 

-1.11i|0-0»l -1.1660-03 2.40'lD-05 1.4450-03 2.5290-02 

2.0850-02 -9-8740-03 -2.0220-02 -1.1810-01 -3.3240+00 
-2.7460-02 1.1510-02 3.7460-02 2.1150-01 1.1100+00 

-6.7210-04 4.4470-03 1.1020-04 -3.7120-03 1.5370-03 

2.0290-06 -2.3790-05 -8.3800-07 1.8700-05 -3.4360-06 


7.0250-05 1.2000-03 -1.2330-03 4.1440-02 5.7810-05 

9.0100-04 1.3860-03 6.8230-04 3.4740-01 -7.9380-04 

3.7140-05 1.9430-06 -8.5780-04 -1.3680-03 -1.7790-05 

-8.9860-03 3.3030-03 1.0260-01 -6.3170-01 -8.3830-03 

6.8650-04 -7.8980-03 2.1970-02 3.9920-01 2.3000-02 

-4.0190-05 -5.8470-05 -1.3210-03 -3.1730-02 -6.5750-04 
4.9870-08 2.3550-07 1.8290-06 4.5770-05 7.2660-07 


In this example cl = 1/E(JA) and c2 = 1/E(JE), where E(JA) and E(JE) denote 
the expected value of the performance indices for the separate flight and 
engine control problems when the previously derived feedback controls are 
used. It should be noted, for this choice of weights, that the gain matrix 
developed for the engine separately is essentially the same as the 
corresponding components of the integrated system feedback matrix. Also, it 
should be noted that each control variable will be active since the gain 
matrix K has all non-zero entries. This structure of the feedback matrix will 
change as the weights in the performance index are varied. The use of a large 
weight on JI negates use of pressure recovery as a control. Variation of the 
parameter c3 would allow one to study the effect of pressure recovery ratio 
upon the integrated system control. 

The introduction of this feedback control improves the stability of the 
overall system and maintains the improvement in short period damping ratio 
which was obtained in the separate airframe control problem. Eigenvalues for 
the open loop integrated system (Int.), for the separate systems with the 
previously derived feedback controls (Air. Ctr., Eng. Ctr.) and for the 
integrated system with feedback control (Int. Ctr.) are listed in tabj.e 6. 
Corresponding short period variables are listed in table 7. 

TABLE 6. Feedback control effects on stability. 

Eigenvalues . 


1 

1 




Airframe 



~i 

1 

_1 

1 

! 

Int 


-2.762E-3 1 

3.743E-4 ± j3.253E-2 | 

-6.782E-1 ± 

j2.201E+0 

1 

1 

Air. 

Ctr. 

-3.242E-3 1 

-6.973E-3 , -1.767E-1 | 

“2.026E+0 + 

jl.258E+0 

1 

1 

1 

Int . 

Ctr. 

-4.729E-3 1 

1 

-1.464E+0 ± jl.ll7E+0 | 

1 

-1.865E+0 ± 

jl .031E+0 

1 

I 

1 

1 




Engine 



1 

_1 

1 

1 

Int 

. 

-5.617E-1 

T j ^ 

1 -1.883E+0 1 -6.585E+0 

, -l.OOOE+1 

1 -1.722E+2 

1 

! 

Eng. 

Ctr. 

-1.484E+0 

1 -3.374E+0 1 -1.007E+1 

± jl.984E+-0 

1 -2.624E+2 

i 

1 

Int. 

Ctr. 

-7.721E-1 

j -3.374E+0 1 -1.007E+1 

± jl.983E+0 

1 -2.624E+2 

1 





J 1 


1 
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TABLE 7. Short period variables. 


! 


I 


cjn - natural frequency (rad/sec) ^ - damping coefficient 

tp - period (sec) th - time to half amplitude (sec) 



wn 

1 

1 5 1 

1 

1 tp 

th 

Int . 

1 

2.303 

1 0.2945 1 

[ 2.729 

1.022 

Air. Ctr. 

2.385 

0.8496 1 

1 2.635 

[ .3422 

Int. Ctr. 

2.130 

1 0.8752 

1 

1 2.948 
1 

.3716 

1 


J 


The transient response of the integrated system to offsets in state 
variables from trim indicates the significance of the above determined 
feedback control. Several general observations can be made concerning the 
integrated system with this feedback control. 

For initial offsets in values of airframe states the following, observations 
can be made. 

1. ) The integrated control system compensates for the coupling between 
airframe and engine variables present in the open loop integrated system 
model and provides good engine control. Recall that small perturbations 
in airframe states can lead to significant changes in engine variables. 

2. ) The closed loop integrated system and the separate closed loop 
airframe model, which was previously analyzed, respond in a similar 
manner. Responses of angle of attack, pitch rate and pitch attitude are 
similar and there appears to be a trade-off between velocity and altitude 
in the two control problems. 

3. ) The dependence of engine states on altitude introduces an implicit 
penalty on deviations of altitude from its trim value. This penalty was 
not present in the separate airframe control problem and is most evident 
in the response of the integrated system to an initial offset in 
altitude. 

4. ) All control variables are active to some extent. 

For initial offsets in values of engine states the following observations can 
be made . 

1. ) Small oscillations in angle of attack, pitch rate and pitch 
attitude, which were present in the open loop system model, are 
amplified. These oscillations are still small relative to the given 
operating point values. 

2. ) The integrated control system provides good engine control. 

3. ) Again all controls are active. 
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The possibility of achieving similar responses with less complex control laws 
will be considered in the next section of this report. 

Transient responses of the open and closed loop integrated model to small 
offsets in the trim values of angle of attack, altitude, compressor speed and 
augmentor pressure are depicted in figures 10 through 13. For an initial 
offset in angle of attack, the system responds with good damping in angle of 
attack, pitch rate and pitch attitude (see figure 10). If this figure is 
compared with figure 5, one observes the similar behavior of angle of attack, 
pitch rate, pitch attitude, glide path angle and thrust. Also, slightly 
better response in velocity is observed for the integrated system. Here the 
horizontal stabilator, 6e, is pulsed for control and it is difficult to see if 
all engine controls are necessary. For an initial offset in altitude, figure 
11, one should notice the increased oscillations of angle of attack, pitch 
rate and pitch attitude and that these disturbances are damped out faster in 
the integrated system than in the separate airframe model (see figure 8). 
Figure 11 shows a rapid altitude change not apparent in figure 8. One should 
note the high h to 6e gain of -5.48 in the integrated system, this gain was 
-0.0037 previously. Other 5e gains are higher, also. The performance index 
for the integrated system has no explicit weight on altitude but altitude 
effects on engine variables are included in the integrated system model. In 
this case the aircraft seems to be reacting to minimize engine perturbations. 
The response of the integrated system to an offset in compressor speed (figure 
12) illustrates coupling present in the integrated sytem with feedback 
control. Similar responses will occur for an offset in fan speed. For 
initial offsets in the trim values of augmentor pressure, fuel flow and 
compressor discharge pressure, very little coupling is apparent in the 
integrated system transient response. This behavior is depicted in figure 13. 

Suboptimal control . Because of the large number of gains present in the 
integrated system feedback matrix, studies were made to determine if 
performance could be maintained using a feedback matrix with a smaller number 
of non-zero gains. Four strategies for determining suboptimal feedback 
matrices were used. The first strategy, Stl, uses a feedback gain matrix 
obtained from the feedback matrices determined in the separate flight and 
engine control problems. In this case the dependence of thrust upon airframe 
states is neglected. In the second strategy, St2, the gain matrix is obtained 
by zeroing gains in the integrated system feedback matrix which correspond to 
control system cross-coupling. In suboptimal strategy three, St3(E), all 
state and control variables, except pitch rate, are normalized relative to 
their given trim value. Gains in the normalized system feedback matrix, whose 
magnitudes are less than some preassigned tolerance s, are set to zero. The 
corresponding gains in the integrated system feedbackmatrix are set to zero. 
In the final strategy, St4(e), sensitivity analysis is used to determine which 
gains should be set to zero. In this method a gain sensitivity matrix 

GSM(g) = (Sen(Xl,g), ..., Sen(X10,g)) 

is calculated, where Xi, i =1, ..., 10 denote the eigenvalues of AAAE - BABE«K 
and g denotes a gain in the feedback matrix K. In this strategy a gain is set 
to zero whenever all entries in the corresponding gain sensitivity matrix are 
less than some preassigned tolerance e. That is, gain g = 0 if 


24 



max{ lSen(Xi,g)| : i = 1,..., 


10 } < £ 


were £ denotes a preassigned tolerance. 

These strategies determine various feedback matrices to be compared with 
the full integrated system feedback matrix previously determined. The 
expected value of the performance index may be used for such a comparison. 
Values of the expected value of performance corresponding to use of different 
feedback matrices are listed in table 8. 

TABLE 8. Suboptimal control. 

Expected performance. 


1 - " I I 



strategy 

E(J) 

1 

1 

T 

1 

Strategy 

E(J) 

i strategy 

E(J) 1 
1 

1 — 

Int.Ctr. 

1 367 

St3( .1) 

1 371 

1 

t St4(.l) 

1 

1 759 t 


Stl 

23 349 

1 

St3( .01) 

1 368 

1 ST4(.01) 

1 376 1 

I 

St2 

3 330 

1 

1 

St3(.001) 

1 367 

( St4(.001) 

1 

1 368 j 
1 


Feedback matrices used in these control strategies are listed in table 9. It 
shoud be noted that gains are being elimated in different orders by strategies 
St3 and St4. 

The performance of the closed loop systems determined by strategies Stl and 
St2 is not as good as that of the system with full feedback control. In each 
case the feedback controls do not compensate, as well as the controls 
determined by the full feedback matrix, for coupling between airframe and 
engine variables. In particular for an offset in the value of an airframe 
state variable, neither of these suboptimal strategies control the engine as 
well as the control determined by the full feedback matrix. A typical 
illustration of this type of behavior is depicted in figure 14, where the 
response of the system with feedback determined by strategy Stl is compared 
with that of the system with full feedback control. 

The number of gains set to zero using strategy St3(£) with e = .1 is 27, 
with £ = .01 is 13 and with e = .001 is 10. In strategy St3 larger values of 
the tolerance lead to unstable feedback systems. The number of gains set to 
zero using strategy St4(£) with e = .1 is 55 , with £ = .01 is 39 and with e = 
.001 is 14. For the feedback control laws determined by strategies St4(.l) 
and St4(.01) the eigenvalues for the ‘closed loop system AAAE - BABE*K compare 
favorably with those of the integrated system with full feedback, these 
eigenvalues are listed in table 10. The corresponding modal matrices for 
these two strategies are listed in table 15 of appendix C. For strategy St4 
larger values of the tolerance yielded eigenvalues which did not compare 
favorably with those of the integrated system with full feedback. 

Since strategy St4(.l) resulted in 55 gains set to zero leaving only 15 
non-zero gains (see table 9), the integrated system with corresponding 
feedback control was analyzed further. In this case, rear compressor variable 
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TABUS 9 . Suboptimal control feedback matrices . 


Strategy Stl 

8 . 113 D -05 

3 . 1530-01 

- 3 . 8870-01 

- 2 . 4780-01 

- 3 . 7180-04 

0.0 

0.0 

0.0 0.0 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9 . 0190-04 

1 . 3950-03 

6 . 5770-04 3 . 4770-01 - 7 . 9350-04 

0.0 

0.0 

0.0 

0.0 

0.0 

3 . 7150-05 

2 . 0320-06 

- 8 . 5800-04 - 1 . 3650-03 - 1 . 7790-05 

0.0 

0.0 

0.0 

0.0 

0.0 

- 8 . 9850-03 

3 . 3290-03 

1 . 0250-01 - 6 . 3130-01 - 8 . 3830-03 

0.0 

0.0 

0.0 

0.0 

0.0 

6 . 7840-04 

- 7 . 9770-03 

2 . 2130-02 3 . 9640-01 2 . 3000-02 

0.0 

0.0 

0.0 

0.0 

0.0 

- 4 . 0200-05 

- 5 . 8620-05 

- 1 . 3210-03 - 3 . 1740-02 - 6 . 5750-04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 0.0 0.0 

Strategy St 2 

4 . 587 D -02 

1 . 7220+00 

- 6 . 4850-01 

- 2 . 7300+00 

- 5 . 4800+00 

0.0 

0.0 

0.0 0.0 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9 . 0100-04 

1 . 3860-03 

6 . 8230-04 3 . 4740-01 - 7 . 9380-04 

0.0 

0.0 

0.0 

0.0 

0.0 

3 . 7140-05 

1 . 9430-06 

- 8 . 5780-04 - 1 . 3680-03 - 1 . 7790-05 

0.0 

0.0 

0.0 

0.0 

0.0 

- 8 . 9880-03 

3 . 3030-03 

1 . 0260-01 - 6 . 3170-01 - 8 . 3830-03 

0.0 

0.0 

0.0 

0.0 

0.0 

6 . 8650-04 

- 7 . 8980-03 

2 . 1970-02 3 . 9920-01 2 . 3000-02 

0.0 

0.0 

0.0 

0.0 

0.0 

- 4 . 0190-05 

- 5 . 8470-05 

- 1 . 3210-03 - 3 . 1730-02 - 6 . 5750-04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 0.0 0.0 


Strategy St 3 (. l ) 


6 e 

4 . 5870-02 

1 . 7220+00 

- 6 . 4850-01 

- 2 . 7300+00 

- 5 . 4800+00 

7 . 0250-05 

1 . 2000-03 

0.0 

0.0 

0.0 

Wfc 

1 . 3640-02 

- 9 . 8220-02 

0.0 

0.0 

5 . 1180-02 

9 . 0100-04 

1 . 3860-03 

0.0 

3 . 4740-01 - 7 . 9380-04 

A 

- 1 . 1140-04 

0.0 

0.0 

0.0 

2 . 5290-02 

3 . 7140-05 

0.0 

- 8 . 5780-04 

0.0 

0.0 

CIW 

2 . 0850-02 

0.0 

0.0 

- 1 . 1810-01 

- 3 . 3240+00 

- 8 . 9880-03 

3 . 3030-03 

1 . 0260-01 

- 6 . 3170-01 - 8 . 3830-03 

RCW 

- 2 . 7460-02 

0.0 

0.0 

2 . 1150-01 

1 . 1100+00 

6 . 8650-04 

- 7 . 8980-03 

2 . 1970-02 

3 . 9920-01 

2 . 3000-02 

BLC 

- 6 . 7210-04 

4 . 4470-03 

0.0 

- 3 . 7120-03 

1 . 5370-03 

- 4 . 0190-05 

- 5 . 8470-05 

- 1 . 3210-03 

- 3 . 1730-02 - 6 . 5750-04 

PR 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


V 

a 

q 

0 

h 

N 1 

N 2 

P 5 

Wf 

P 2 



TABLE 9 


concluded . 


Strategy St3 (. 01 ) 

4 . 587D-02 1 . 722D+00 

- 6 . 4850-01 

- 2 . 7300+00 

- 5 . 4800+00 

7 . 0250-05 

1 . 2000-03 

- 1 . 2330-03 

4 . 1440-02 

5 . 7810-05 

1 . 364D-02 - 9 . 822D-02 

0.0 

5 . 6800-02 

5 . 1180-02 

9 . 0100-04 

1 . 3860-03 

6 . 8230-04 

3 . 4740-01 

- 7 . 9380-04 

- 1 . 114D-04 - 1 . 166D-03 

0.0 

1 . 4450-03 

2 , 5290-02 

3 . 7140-05 

1 . 9430-06 

- 8 . 5780-04 

- 1 . 3680-03 

- 1 . 7790-05 

2 . 085D-02 " 9 . 874D-03 

- 2 . 0220-02 

- 1 . 1810-01 

- 3 . 3240+00 

- 8 . 9880-03 

3 . 3030-03 

1 . 0260-01 

- 6 . 3170-01 

- 8 . 3830-03 

- 2 . 746D-02 1 . 151D-02 

0.0 

2 . 1150-01 

1 . 1100+00 

6 . 8650-04 

- 7 . 8980-03 

2 . 1970-02 

3 . 9920-01 

2 . 3000-02 

- 6 , 7210-04 4 . 4470-03 

1 . 1020-04 

- 3 . 7120-03 

1 . 5370-03 

- 4 . 0190-05 

- 5 . 8470-05 

- 1 . 3210-03 

- 3 . 1730-02 

- 6 . 5750-04 

0.0 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


Strategy St4 (. l ) 


4 . 587D-02 

1 . 7220+00 

- 6 . 4850-01 

- 2 . 7300+00 

- 5 . 4800+00 

0.0 

1 . 2000-03 

0.0 

0.0 

0.0 

1 . 3640-02 

- 9 . 8220-02 

0.0 

5 . 6800-02 

0.0 

9 . 0100-04 

1 . 3860-03 

0.0 

3 . 474D-01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 . 7140-05 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 8 . 9880-03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 6 . 575D-04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


Strategy St4 {. 01 ) 


6e 

4 . 5870-02 

1 . 7220+00 

- 6 . 4850-01 

- 2 . 7300+00 

- 5 . 4800+00 

0.0 

1 . 2000-03 

0.0 

4 . 1440-02 

0.0 

Wfc 

1 . 3640-02 

- 9 . 8220-02 

- 6 . 1780-03 

5 . 6800-02 

0.0 

9 . 0100-04 

1 . 3860-03 

0.0 

3 . 4740-01 

0.0 

A 

- 1 . 1140-04 

- 1 . 1660-03 

0.0 

1 , 4450-03 

2 . 5290-02 

3 . 7140-05 

0.0 

- 8 . 5780-04 

0.0 

0.0 

CIW 

0.0 

0.0 

0.0 

0.0 

0.0 

- 8 . 9880-03 

0.0 

0.0 

0.0 

0.0 

RCW 

- 2 . 7460-02 

0.0 

0.0 

0.0 

0.0 

0.0 

- 7 . 8980-03 

0.0 

0.0 

0.0 

BLC 

- 6 . 7210-04 

4 . 4470-03 

0.0 

- 3 . 7120-03 

0.0 

- 4 . 0190-05 

- 5 . 8470-05 

- 1 . 3210-03 

- 3 . 1730-02 

- 6 . 5750-04 

PR 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 


V 

a 

q 

e 

h 

N1 

N2 

P5 

Wf 

P2 


lo 



vane position, RCW, and pressure recovery ratio, Pr, are not used as controls 
and there is only one gain associated with nozzle area, a, inlet guidevane, 
CIW, and compressor bleed, BLC. Also changes in compressor speed are fed 
back to control 6e, whereas, changes in velocity, angle of attack and pitch 
attitude are fed back to control fuel flow. Transient responses of the 
integrated system using suboptimal feedback matrices determined by strategy 
St4(.l) and the full feedback matrix previously determined are nearly the same 
for offsets in any state variable. Typical responses for the integrated 
system using the control determined by strategy St4(.l) and the full feedback 
matrix previously determined are depicted in figures 15 and 16. In each case, 
it should be observed that the integrated system with suboptimal feedback 
behaves in a manner very similar to the system with full feedback control. 

TABLE 10. Suboptimal control effects on stabilty. 

System eigenvalues. 


j Airframe I 

1 1 

1 

Int. Ctr. 

-4.729E-3 

-1.464E+0 ± jl.ll7E+0 

1 -1.865E+0 

± 

jl.031E+0 

n 

1 

i 

ST4(.l) 

-4.800E-3 

-1.395E+0 ± jl.l05E+0 

j -1.893E+0 

+ 

jl.083E+0 

1 

u 

ST4(.01) 

-4.735E-3 



-1.450E+0 ± jl.lOOE+0 

[ -1.863E+0 

.. 1 

+ 

jl.052E+0 

1 

1 

[ Engine j 

1 _ _ [ 

1 

Int. Ctr,. 

-7.721E-1 

1 

-3.374E+0 1 -1.007E+1 

± jl.983E+0 

1 

-2.624E+2 

1 

i 


ST4( .1) 

-7 .105E-1 

-3.519E+0 1 -0.902E+1 

± j0.803E+0 

1 

-2.590E+2 

1 

1 

L_ 

ST4( .01) 

-7.689E-1 



“3.352E+0 1 -1.015E+1 

1 

± jl.982E+0 

1 

1 

-2.577E+2 

1 

J 
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CONCLUDING REMARKS 


In this report numerical and analytical techniques are used to develop 
models of aircraft subsystem interactions and to analyze the effect of 
subsystem interactions upon system stability and control. 

The inclusion of subsystem interactions in the airframe/propulsion system 
model changed the overall system behavior. Also, the approach used to 
develop the flight/propulsion system model in this report is general 
since standard flight and propulsion system models are combined to form 
the integrated system model. 

Methods presented in this report provide a means to identify interaction 
parameters in an operating point model of a flight/propulsion model. 

The model, developed in this report, provides a means for further 
parameter studies. For example, the effect of changing the location of 
the engine below the center of gravity or the effect of using inlet 
pressure recovery ratio as a control parameter could be readily analyzed 
using this model. 

The problem of developing adequate models of subsystem interactions 
remains difficult. For example, the effect of variable geometry inlets 
upon aerodynamic forces and moments have not been included in the model . 
Also angle of attack and sideslip angle effects upon engine airflow have 
not been considered. 

Linear quadratic regulator methods and numerical linear algebra 
techniques have provided flexible means to analyze the control of 
integrated systems. 

For the separate airframe model, feedback control provides substantial 
improvement in short period damping. 

For the integrated flight/propulsion model, feedback control conpensates 
for coupling present in the model and provides good overall system 
stability. In the model example, the aircraft appears to be reacting to 
minimize engine perturbations. This behavior results from the dependence 
of engine variables upon Mach number and altitude. In the performance 
index used for this example, the engine was referenced to a fixed trim 
point. Perhaps the engine should have been referenced to trim values 
which change with Mach number and altitude. 

Analysis of suboptimal control strategies, indicates that performance of 
the closed loop integrated system can be maintained with a feedback 
matrix in which the number of nonzero gains is small relative to the 
number of components in the feedback matrix. 

A method based on senitivity analysis, proved to be an effective means 
for determining which gains in the integrated system feedback matrix can 
be set to zero while, at the same time, maintaining system performance. 
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APPENDIX A 


AIRCRAFT EFFECTS ON ENGINE 


The off-design performance model for a turbofan engine developed by F. J. 
Lallman [21] expresses total temperatures and pressures throughout the engine 
along with specific thrust and fuel to air ratio as functions of flight 
condition (Mach number and ambient temperature), inlet performance (pressure 
recovery ratio), fan pressure ratio, compressor pressure ratio and burner 
temperature. This model may be used to obtain engine perturbation equations 
which describe the effect of varying flight condition and inlet performance 
upon propulsion system operation. 

Notation consistent with that of reference [21] will be used in this 
appendix. For a given operating point, perturbation equations will be 
presented which express the effect of changes in Mach number,^ altitude amd 
pressure recovery ratio upon engine temperatures, pressures and performance 
variables . 


The engine characteristics are as listed in Table I of [21] . A full power 
engine operating point is assumed, this condition corresponds to a fan 
pressure ratio of 2.9, a compressor pressure ratio of 7.93 and a burner 
temperature of 1559 K. The engine design point model is as follows; 


6 


“s 


/ 5.256 



0 < ho < 11 


( 0.2233 


^-0.157(ho-ll) 


11 < ho < 20 


60 = O./0Os)^'^^® 6.S 

61 - ^10^0 

/ 

= 2.9 

62 = 7.93 61, 

63 = 0.95 62 

4 290 

6^ = (1.016 - 0.3908 0i) 63 


Sb = 6. 

and 


00 


s 


1 -0.02256ho 
0.7519 


0 < ho < 11 
11 < ho < 20 
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00 = 

(1 + O.199m/)0„^ 


II 

HI 

o 

00 


CD 

II 

1.38701 


&2 

1.86401, 


03 = 

5.41 


04 = 

5.488 - 1.90301 


05 = 

4.677 - 1.23401 


Vo 

34O.4(0Og)'/^M„ 


WA 

79.8876i/(0i)^^^ 


WF/WA = 

0.02332 - 0.0108601 


WF 

(WF/WA )WA 


ST 

STl + ST2 - Vp 


STl 

482.9(l-(6s/6og)‘°‘^^^® 

1/2 

) ‘ (0s) 

ST2 

273.9(l-(6i,/6og)"°'^^^ 

1/2 

) ^ (01. 

TH 

ST»WA 



where the 6's represent engine pressures relative to sea level standard (101.4 
kPa) , the 0's represent engine temperatures relative to sea level standard 
(288.2 K) , WA fan airflow, WF/WA fuel to air ratio, v“ aircraft velocity, ST 
specific thrust (thrust per unit mass flow), and TH denotes thrust. 

Using the notation 


Aho=h, Ar^ Q=Pr, AMc,=M, AWF=Wf ATH=Th 


the perturbation equations for this design point are as follows: 
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0 < ho < 11 


and 


and 


AVp 

AWA 


4 256 

I 5.2560og A0og 


A6„^ 


I -0.15776ogh 

A6o 

= 

Oa/0Os)^'^^®A«.g + 3.5186 

A6i 

= 

qA5q ■1'6oP^ 

/ 

A6i, 

= 

2.9A6i 

A62 

= 

7.93A6i, 

A63 

= 

O. 95 A 62 

A64 

= 

-1.67656 3(1.016-0.390881)^ 

A65 

= 

A64 


11 < ho < 20 


2.518 


0/ ^D, 


4.29 


A6. 




( -0.02256h 0 

< 

ho 

< 

11 

A0og 


(0 11 

< 

ho 

< 

20 

A0O 

=r 

O.3980o^Mo»M + (1 + 

0, 

,199Mo 

2 

) 

A01 

= 

A0O 





A01 , 

= 

1 .387A01 





A02 

= 

1.864A01 , 





A03 

= 

0.0 





.9- 

CD 

< 

= 

-1.9O3A01 





A05 

= 

-1 .234A01 






17O.2Mo(0Og)''^^A0Og + 34O.4(0Og)'^^M 


= WA(A6i/6i - O.5A0i/0i) 
A(WF/WA) = -O.O1O86A01 

Wf = WAA(WF/WA) + (WF/WA)AWA 
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AST = ASTI + AST2 - Av^ 

ASTI = 65.14305-®((l-(55/6.g)"'^^^®>'’^ 


+ 241.4505"'^ (1 -(65/fiog)‘‘^^^®)’^A05 

AST2 = 36.429(01, )’^(1-(Si,/6o 5)’’^^^)"'^ ( 6i , /6 '^®®A(5 i , /6,^) 

+ 136. 9501, ■•^(l-<6i,/6,^)‘‘^®^)‘^ A01,. 


Th = (AST)»WA + (AWA)*ST 

For the control design example under consideration the nominal values of 
Mach number, altitude and inlet pressure recovery ratio are 


Mj, = 0.9 ho = 12. km rit, = 1.0 

f 

The engine variables of interest are P2 , P5 , Wf, Th and Wa where 

P2 = 101.462 and P5= 101.46s. 

For this flight condition and pressure recovery ratio the resulting 
perturbation equations may be written as 


P5 = 

42.02 

M - 

15.938 

h 

+ 

101.016 

Pr 

Wf = 

2.717 

M - 

0.5905 

h 

+ 

3.7487 

Pr 

P2 = 

635.55 

M - 

90.631 

h 

+ 

554.812 

Pr 

Th = 

8267.0 

M - 

2121.0 

h 

+ 

17826.0 

Pr 

Wa = 

25.678 

M “ 

4.350 

h 

+ 

27.586 

Pr 


These equations describe perturbations in P5, Wf, P2, Th, and Wa due to 
perturbations in Mach number, altitude, and pressure recovery ratio necessary 
to maintain the operating condition. 

The engine dynamics about the operating point are described by the 
differential equations 

Dxe/Dt = AE*xe + BE»ue + CEA»ya + CEI*yi 
ye = FE*xe + GE*ue + HEA»ya + HEI*yi. 
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In order to determine the interaction matrices CEA, CEI, HEA, HEI, it is 
convenient to rewrite the above perturbation equations in matrix form as 
follows 




M 


M 

n tviT 

M 

= El* 

h 

[Wf] = E2* 

h 

1 W;:) 1 

h 


Pr 


Pr 


Pr 


where El, E2, and E3 are appropriate matrices. 


At the given operating point it is assumed that ya and yi have no direct 
effect upon the variables N1 , N2, and Wf. Consequently the first, second, and 
fourth rows of CEA and CEI are zero v/hile the remaining entries of CEA, CEI, 
HEA, and HEI are chosen so that the steady state response of the dynamic 
system to changes in ya and yi agrees with the previously derived perturbation 
equations . 


Algebraically, these assumptions determine the following equations for the 
unknown entries of the matrices CEA, CEI, HEA, and HEI. With no control the 
steady state response of the engine is given by the equations 


xe = -(AE)"^[ CEA-ya + CEI»yi ] 
ye = FE*xe + HEA«ya + HEI*yi. 


In order that these equations agree with the previously determined 
perturbation equations for P5, P2, Th and Wa it is necessary that 


El 


0 . 0 . 0 . 1 . 0 . 0 . 

0 . 0 . 0 . 0 . 0 . 1 . 


•AE 


-1 


[CEA, CEI] 


•AE"^[CEA, CEI] 


v/here [CEA, CEI] and [HEA, HEI] denote the matrices whose columns are the 
corresponding columns of CEA, CEI, HEA, and HEI. FEl and FE2 denote the first 
and second rows, respectively, of the matrix FE. Using the modelling 
assumption that the first two rows of CEA and CEI are zero, these equations 
may be inverted to obtain the unknown components of CEA, CEI, HEA, and HEI. 


[HEA, HEI] = E3 + 


FEl 


FE2 
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APPENDIX B 


EIGENVALUE SENSITIVITY CALCULATIONS 

In the case of distinct eigenvalues the real Jordan canonical form A of a 
matrix A is a block diagonal matrix. The i-th diagonal block Ai of A is if 
X is a real eigenvalue of A or 



if X = a + jb is a complex eigenvalue. The modal matrix T such that 

A«T = T»A (B.l) 

is determined by the eigenvectors of A. If X is a real eigenvalue with 

corresponding eigenvector u then the corresponding column of T equals u. 

Whereas if X is complex with eigenvector u + jv then the corresponding columns 
of T equal u and v. In the real case the fundamental mode corresponding to X 
is 

X(t) = u*e^^. 

In the complex case the fundamental modes corresponding to X = a + jb are 

X(t) = u*cos(bt) - v«sin(bt) ] 

Y(t) = u«sin(bt) + v«cos(bt) ]. 

If the matrix A depends on a parameter s then differentiation of equation 
(B.l) yields 

T”^*DA/De*T + A*T**^»DT/D£ = t“^*DT/D£*A + DA/De . 

In this last equation if the i-th diagonal block of A corresponds to a real 
eigenvalue then 

( A*t'^»DT/De ). . = ( t"^*DT/De«A ). . 

J. , 1 1 ± 

and it follows that 

DX/De = ( t“^*DA/De»T ). .. 

If the i-th diagonal block of A corresponds to a complex eigenvalue then 
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(A*t'^*DT/De). . 

± f ± 

= 

(a*T~^*DT/DE)^ ^ -t (b*T“^*DT/D£)^_j^^^^ 

(A.T-^-DT/Dc>._i^^ 


(a*T“^»DT/DE)^^^^^ + (b»T"^*DT/DE)jL+i i+1 

(A*T’^.DT/De).^j^^^ 

= 

(-b»T"^*DT/DE)^ ^ + (a*T“^«DT/De)^^^^^ 

(A.T-^.DT/De).^1 



(t'^«DT/De*A) . . 

J. f X 

= 

(a*T“^*DT/D£) . . - (b*T‘^»DT/De) . . 

X ^ X X , i * j. 

(T-^DT/De.A)i 

= 

(b*T~^*DT/D£)^ ^ + (a»T"^*DT/D£)^ 

(T-^DT/De.A).^l_i 

= 

(a.T'^.DT/Dc),^^ . - (fT'^*DT/D£).^j 

(T-^.DT/DE.A)i^l_..,l 

= 


Since 




t'^^DA/De^T + A»T ^•DT/De = T ^•DT/De*A + DA/De 
t follows that 

(t"^*DA/De«T)^ ^ + (t“^*DA/De*T)^^^ + (A»T ^•DT/De)^^^ 

= (t'^»DT/De»A)^ ^ + 2Da/Dc 

and consequently 

Da/DE = ((t"^*DA/De*T)^ ^ + <t'^*DA/De*T)^_^^^^^^)/2 

Similarly 

Db/Ds = ((t"^»DA/De*T)^ - (t“^*DA/DeT)^^^^^)/2. 

These identities .are used to calculate sensitivity derivatives in this report. 
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APPENDIX C 


SYSTEM MATRICES 

System matrices for the models presented in this report are contained in 
this appendix. In table 11 matrices associated with the separate airframe 
model are listed. This data was supplied by NASA. In table 12 matrices 
associated with the separate engine model are listed. This data was derived 
from data presented in the report [15] . System matrices for the integrated 
flight/propulsion model are presented in table 13. Modal matrices for the 
separate flight and engine models and for the integrated model, are listed in 
table 14. In table 15 the eigenvalues and modal matrices for suboptimal 
control strategies Stl, St2, St3 and St4 are listed. The modal matrices 
listed in tables 14, 15 have been denoted by the symbol T in the body of this 
report. The colunms of these matrices consist of eigenvectors corresponding 
to real eigenvalues and the real and imaginary part of complex eigenvectors 
corresponding to complex eigenvalues. The order in which real and imaginary 
parts of eigenvectors appear in these modal matrices is the order in which 
eigenvalues are listed in tables 14, 15. 
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tJ 

cx> 


— 


AA 


— 

- 1 . 57290-02 

- 1 . 19230+01 

0.0 

- 9 . 80650+00 

2 . 00290-01 

- 3 . 67560-04 

- 5 . 97400-01 

1 . 00000+00 

0.0 

5 . 86080-03 

- 3 . 17330-03 

- 4 . 84240+00 

- 7 . 45220-01 

0.0 

5 . 04860-03 

0.0 

0.0 

1 . 00000+00 

0.0 

0.0 

0.0 

- 2 . 65560-01 

0.0 

2 . 65560-01 

0.0 


- 

1 : 22460-04 

0.0 

CAE 

0.0 

0.0 

0.0 

0.0 


- 3 . 51750-08 

0.0 

0.0 

0.0 

0.0 

0.0 


2 . 55610-06 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


3 . 38870-03 

0.0 

0.0 

0.0 

FA 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 00000+00 



1 . OOOOD +00 

0.0 

FFA 

0.0 

0.0 

0.0 

0.0 

- 1 . 00000+00 

0.0 

1 . 00000+00 

0.0 

0.0 

0.0 

1 . 00000+00 

0.0 

0.0 


4 . 00000-08 

0.0 

0.0 

WA 

0.0 

2 . 50000-01 

0.0 

0.0 

0.0 

1 . 00000+00 



BA -1 
- 2 . 3459 D +00 
- 6 . 4267 D -02 
- 7 . 4235 D +00 
0.0 
0.0 


GGA 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


- UA 

4 . 0000 D +00 0.0 

0.0 2 . 5000 D 


-n_ 


TABLE 11. Airframe matrices 



-9.606D-01 

-5.9660-01 

AE 

-1.6110+02 

1 .0820+03 

1.8810+01 

7.955D-01 

-1 .6440+00 

-3.8410+01 

-1.5460+02 

9.8230+00 

1.786D-02 

-3.5720-02 

-8.8860+00 

4.1260+01 

5.7560-01 

0.0 

0.0 

0.0 

-1.0000+01 

0.0 

j;;5.128D+00 

1.2520+01 

6.5720+02 

9. 6260+03 

-1.6970+Q2i 

-1. 9360+02 

-3.8860+03 

BE 

-2.7230+01 

-8.0010+00 

2.8000+03 

3. 3290+01 

9.1360+02 

1.7470+01 

-1.4170+01 

7.1940+03 

-3.6540+00 

-1.541 0+03 

6.2880-01 

-5.1300-01 

-2. 5590+01 

1.0000+01 

0.0 

0.0 

0.0 

0.0 

jJ. 6400+03 

-1.4370+04 

-8.5430+01 

9.7290+01 

-1.3270+0^, 

1 .6490+00 

-2. 0280+00 

FE 

-2.7850-02 

-2.7050+02 

2.7790+01 

1 .9620-03 

8.6410-07 

-8.6710-04 

8.1530-02 

2.8780-07 

8.3060-02 

-5.0020-02 

-2.9120+00 

1.8870+02 

-8.5330-02 

2.0680-04 

1.8620-06 

-8.7330-03 

9.4280-03 

4.2630-05 

-6.2040-05 

1 .5240-04 

6.9600-03 

6.7220-02 

-1.6060-03 

5.7050-05 

-3.7290-05 

-6.3700-03 

-3.5380-02 

1.1570-OJi. 

1 .8820+01 

1.1410+04 

GE 

4.0250+01 

-1.3780+01 

1.5030+04 

2.3330-04 

9.9600-01 

8.1600-02 

1.6980-04 

7.0130-01 

6. 2000+00 

3.0670+02 

1.3570+00 

-2.4410-01 

1 . 0720+03 

7.7750-04 

2. 7700-01 

7.3550-04 

1.4060-04 

5.9090-01 

-2.0630-02 

-2. 1430-01 

-1.0700-03 

-1.6210-03 

-4.0190-01 

1.1480-02 

3.0740-01 

1.9010-03 

-2.3970-04 

2.6330-01 



■ — 

CEA 




CEI 


HEA 




HEl 




0.0 


0.0 



0.0 

-5.5570+03 -4.6070+02 



7.2930+03 



0.0 


0.0 



0.0 

2.0210+01 -5.8780+00 



3.7260+01 



1.0660+02 - 

9.5940+01 


1.0660+02 

0.0 

0.0 




0.0 




0.0 


0.0 



0.0 

0.0 

0.0 




0.0 




4..2370+04 -3.5050+03 


4.2370+04 

0.0 

0.0 




0.0 










0.0 

0.0 




JO.O 

— 


— 




WE 






UE 




— 

5.0000-04 

0.0 


0.0 

0.0 

0.0 

6.5500+01 

0.0 

0 

.0 

0.( 

D 

0 

0 

0.0 


5.0000-04 

0.0 

0.0 

0.0 

0.0 

4.6340+04 

0 

.0 

0.0 

0 

0 

0.0 


0.0 


1.0520-02 

0.0 

0.0 

0.0 

0.0 

5 

.0000-01 

0.( 

1 

0 

0 

0.0 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2.0000-01 

0 

0 

0.0 


0.0 


0.0 

0.0 

2.1040-02 

0.0 

0.0 

0 

.0 

0.0 

1 

0000+Qi 


OJ 

VO 


TABLE 12. Engine system matrices 



AAAE 


- 1 . 8040-02 

- 1 . 1920+01 

0.0 

- 9 . 8060+00 

1 . 4390-01 

2 . 0190-04 

- 2 . 4840-04 

- 3 . 4110-06 

- 3 . 3130-02 

3 . 4040-03 

- 3 . 669 D -04 

- 5 . 9740-01 

1 . 0000+00 

0.0 

5 . 8770-03 

- 5 . 7990-08 

7 . 1330-08 

9 . 7970-10 

9 . 5160-06 

- 9 . 7760-07 

- 3 . 2210-03 

- 4 . 8420+00 

- 7 . 4520-01 

0.0 

3 . 8710-03 

4 . 2140-06 

- 5 . 1840-06 

- 7 . 1190-08 

- 6 . 9150-04 

7 . 1040-05 

0.0 

0.0 

1 . 0000+00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 2 . 6560-01 

0.0 

2 . 656 D -01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 9 . 6060-01 

- 5 . 9660-01 

- 1 . 6110+02 

1 . 0820+03 

1 . 8810+01 

0 . 0 , 

0.0 

0.0 

0.0 

0.0 

7 . 9550-01 

- 1 . 6440+00 

- 3 . 8410+01 

- 1 . 5460+02 

9 . 8230+00 

3 . 6110-01 

0.0 

0.0 

0.0 

- 9 . 5940+01 

1 . 7860-02 

- 3 . 5720-02 

- 8 . 8860+00 

4 . 1260+01 

5 . 7560-01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 1 . 0000+01 

0.0 

1 . 4360+02 

0.0 

0.0 

0.0 

- 3 . 5050+03 

- 5 . 1280+00 

1 . 2520+01 

6 . 5720+02 

9 . 6260+03 

- 1 . 6970+02 


BABE 


- 2 . 3460+00 

2 . 3050-03 

1 . 3970+00 

4 . 9290-03 

- 1 . 6870-03 

1 . 8410+00 

8 . 9320-01 

- 6 . 4270-02 

- 6 . 6210-07 

- 4 . 0120-04 

- 1 . 4160-06 

4 . 8460-07 

- 5 . 2870-04 

- 2 . 5650-04 

- 7 . 4240+00 

4 . 8120-05 

2 . 9150-02 

1 . 0290-04 

- 3 . 5210-05 

3 . 8420-02 

1 . 8640-02 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 1 . 9360+02 

- 3 . 8860+03 

- 2 . 7230+01 

- 8 . 0010+00 

2 . 8000+03 

0.0 

0.0 

3 . 3290+01 

9 . 1360+02 

1 . 7470+01 

- 1 . 4170+01 

7 . 1940+03 

0.0 

0.0 

- 3 . 6540+00 

- 1 . 5410+03 

6 . 2880-01 

- 5 . 1300-01 

- 2 . 5590+01 

1 . 0660+02 

0.0 

1 . 0000+01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 1 . 6400+03 

- 1 . 4370+04 

- 8 . 5430+01 

9 . 7290+01 

- 1 . 3270+05 

4 . 2370+04 


FAFE 


3 . 3890-03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 0000+00 

0.0 

O'.O 

0.0 

0.0 

0.0 

1 . 8830+01 

0.0 

0.0 

0.0 

- 4 . 6070+02 

1 . 6490+00 

- 2 . 0280+00 - 2 . 7850-02 

- 2 . 7050+02 

2 . 7790+01 

6 . 8480-02 

0.0 

0.0 

0.0 

- 5 . 8780+00 

1 . 9620-03 

8 . 6410-07 - 8 . 6710-04 

8 . 1530-02 

2 . 8780-07 

0.0 

0.0 

0.0 

0.0 

0.0 

8 . 3060-02 

- 5 . 0020-02 - 2 . 9120+00 

1 . 8870+02 

- 8 . 5330-02 

0.0 

0.0 

0.0 

0.0 

0.0 

2 . 0680-04 

1 . 8620-06 - 8 . 7330-03 

9 . 4280-03 

4 . 2630-05 

0.0 

0.0 

0.0 

0.0 

0.0 

- 6 . 2040-05 

1 . 5240-04 

6 . 9600-03 

6 . 7220-02 

- 1 . 6060-03 

0.0 

0.0 

0.0 

0.0 

0.0 

5 . 7050-05 

- 3 . 7290-05 - 6 . 3700-03 

- 3 . 5380-02 

1 . 1570-04 


GAGE 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 8820+01 

1 . 1410+04 

4 . 0250+01 

- 1 . 3780+01 

1 . 5030+04 

7 . 2930+03 

0.0 

2 . 3330-04 

9 . 9600-01 

8 . 1600-02 

1 . 6980-04 

7 . 0130-01 

3 . 7260+01 

0.0 

6 . 2000+00 

3 . 0670+02 

1 . 3570+00 

- 2 . 4410-01 

1 . 0720+03 

0.0 

0.0 

7 . 7750-04 

2 . 7700-01 

7 . 3550-04 

1 . 4060-04 

5 . 9090-01 

0.0 

0.0 

- 2 . 0630-02 

- 2 . 1430-01 

- 1 . 0700-03 

- 1 . 6210-03 

- 4 . 0190-01 

0.0 

0.0 

1 . 1480-02 

3 . 0740-01 

1 . 9010-03 

- 2 . 3970-04 

2 . 6330-01 

0.0 


TABLE 13. Integrated system matrices. 



FF 


1 . OOOD +00 

0.0 

0.0 

0.0 

0.0 

0.0 

- 1 . 0000+00 

0.0 

1 . 0000+00 

0.0 

0.0 

0.0 

1 . 0000+00 

0.0 

0.0 

1 . 8830+01 

0.0 

0.0 

0.0 

- 4 . 6060+02 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 6490+00 

- 2 . 0280+00 

- 2 . 7850-02 

- 2 . 7050+02 

2 . 7790+01 

1 . 0000+00 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 0000+00 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 0000+00 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 0000+00 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 0000+00 






GG 






0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

1 . 8820+01 

1 . 1410+04 

4 . 0250+01 

- 1 . 3780+01 

1 . 5030+04 

7 . 2930+03 



0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


f— 




w 




— 

1 . 2740-08 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 . 9620-02 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 . 1850-01 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

7 . 9620-12 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 6890-04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 6890-04 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 . 5530-03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 . 1070-03 


1 . 2740+00 

0.0 

0.0 

U 

0.0 

0.0 

0.0 

0.0 



0.0 

2 . 2130+01 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

1 . 5660+04 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

1 . 6890-01 

0.0 

0.0 

0.0 



0.0 

0.0 

0.0 

0.0 

6 . 7570-02 

0,0 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

3 . 3780+03 

0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 000 D +06 



TABLE 13. concluded 



to 


Eigenvalues 


Airframe 


Engine 


Re. 

- 2 . 762 D -03 
3 . 743 D -04 
3 . 743 D -04 
- 6 . 782 D -01 
- 6 . 782 D -01 
- 5 . 617 D -01 
- 1 . 8840+00 
-6. 5850+00 
-1.0000+01 
- 1 . 7220+02 


Im. 

0.0 

3 . 2530-02 

- 3 . 2530-02 

2 , 2000+00 

-2.2000+00 

0.0 

0.0 

0.0 

0.0 

0.0 


Separate systems modal matrix. 


1 . 0000+00 

- 9 . 9960-01 

0.0 

8 . 6880-01 

- 4 . 1360-01 

0.0 

0.0 

0.0 

0.0 

0.0 

6 . 6030-04 

6 . 9910-04 

- 1 . 1240-05 

- 3 . 7610-03 

- 1 . 0450-01 

0.0 

0.0 

0.0 

0.0 

0.0 

1 . 6930-06 

- 1 . 0460-04 

4 . 4570-05 

2 . 3070-01 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6 . 1280-04 

1 . 3330-03 

3 . 2300-03 

- 2 . 9500-02 

- 9 . 5730-02 

0.0 

0.0 

0.0 

0.0 

0.0 

4 . 5730-03 

2 . 6520-02 

- 4 . 8710-03 

1 . 8460-03 

2 . 5380-03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 . 8050-01 

- 8 . 8930-01 

9 . 7550-01 

9 . 9770-01 

- 1 . 1250-01 

0.0 

0.0 

0.0 

0.0 

0.0 

8 . 1290-01 

4 . 5290-01 

- 9 . 5010-02 

3 . 7560-02 

- 5 . 7410-02 

0.0 

0.0 

0.0 

0.0 

0.0 

9 . 5850-04 

6 . 3890-04 

5 . 6590-02 

4 . 5630-02 

- 3 . 4970-03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 2 . 0890-03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

4 . 6300-02 

6 . 3470-02 

1 . 9000-01 

3 . 2790-02 

9 . 9200-01 


Eigenvalues 


Airframe 


Engine 


Re. 

/ 1 . 9120-03 
1 - 3 . 6540-04 
< - 3 . 6540-04 
I - 6 . 7810-01 
^ - 6 , 7810-01 
- 5 . 6280-01 
- 1 . 8830+00 
- 6 . 5870+00 
-1. 0000+01 
- 1 . 7220+02 


Im. 

0.0 

3 . 6470-02 

- 3 . 6470-02 

2.2000+00 

-2.2000+00 

0.0 

0.0 

0.0 

0.0 

0.0 


Integrated system modal matrix. 


6 . 3110-02 
- 4 . 0400-05 
- 7 . 9290-08 
- 4 . 1460-05 
- 1 . 4750-04 
4 . 8440-01 
8 . 5960-01 
1 . 1330-02 
0.0 

1 . 4910-01 


- 6 . 5670-02 

4 . 3300-05 

- 8 . 8240-06 

4 . 0790-05 

1 . 7910-03 

7 . 6140-01 

- 5 . 5720-01 

- 3 . 7880-02 

0.0 

- 3 . 0340-01 


6 . 5430-03 
- 4 . 3370-06 
1 . 3990-06 
2 . 4160-04 
2 . 9240-07 
3 . 6970-02 
1 . 0150-01 
1 . 1470-03 
0.0 

1 . 6410-02 


1 . 5780-01 

- 5 . 9290-04 

4 . 1860-02 

- 5 . 2730-03 

3 . 3290-04 

2 . 4170-01 

- 3 . 8010-01 

7 . 4520-03 

0.0 

1 . 1810-01 


- 7 . 4330-02 

- 1 . 8980-02 

1 . 9710-04 

- 1 . 7400-02 

4 . 6220-04 

- 5 . 3620-01 

- 6 . 5150-01 

- 2 . 2620-02 

0.0 

- 1 . 9420-01 


- 1 . 2370-04 
3 . 9450-07 
- 1 . 2010-08 
2 . 1340-08 
1 . 7610-07 
5 . 8080-01 
8 . 1280-01 
9 . 3810-04 
0.0 

4 . 6090-02 


3 . 7940-05 
- 2 . 7170-07 
3 . 4160-07 
- 1 . 8140-07 
- 1 . 2730-08 
- 8 , 8940-01 
4 . 5270-01 
6 . 4570-04 
0.0 

6 . 3520-02 


- 1 . 3040-04 

4 . 9420-07 

- 2 . 7590-06 

4 . 1890-07 

3 . 0380-09 

9 . 7560-01 

- 9 . 4790-02 

5 . 6590-02 

0.0 

1 . 8990-01 


- 3 . 7150-05 

9 . 5400-08 

- 8 . 0360-07 

8 . 0360-08 

3 . 9950-10 

9 . 9770-01 

3 . 7580-02 

4 . 5630-02 

- 2 . 0880-03 

3 . 2790-02 


- 1 . 9560-05 
7 . 9880-09 
- 4 . 1020-07 
2 . 3830-09 
8 . 6460-12 
-1 .1250-01 
- 5 . 7410-02 
- 3 . 4970-03 
0.0 

9 . 9200-01 


TABLE 14 . Modal matrices 



Stl eigenvalues 

i 


Airframe 


Engine 


Re. 

-3.076D-03 
-7.380D-03 
-1 .744D-01 
-2.019D+00 
-2.019D+00 
. -1.483D+00 
I -3.374D+00 
< -1 .007D+01 
I -1.007D+01 
' -2.624D+02 


Im. 

0.0 

0.0 

0.0 

1 . 268D+00 
-1.2680+00 
0.0 
0.0 

1 .9840+00 
-1.9840+00 
0.0 


Stl modal matrix 


1.6120-01 

-9.4250-05 

3.0440-07 

-9.8980-05 

4.0850-04 

2.1040-01 

9.6080-01 

3.4210-04 

-1.0880-03 

8.1610-02 


-1.2790-02 
3.8260-05 
2.8320-07 
-3. 8380-05 
2.7580-03 
9.6900-01 
-1.9660-01 
-2.7430-02 
-5.1800-04 
-1.4620-01 


-2.9530-02 

5.6030-04 

1.8780-04 

-1.0770-03 

2.4930-03 

9.3900-01 

-3.0820-01 

-2.5400-02 

-3.8890-04 

-1.4780-01 


- 5 . 5280-02 
-1.3280-02 
2.9700-02 
-1 .1820-02 
5.7540-05 
-2.3490-01 
-4.8820-01 
-3.9050-03 
7.0440-04 
-2.9730-02 


-1.0530-01 

-7.8920-03 

-5.7130-03 

-4.5940-03 

-3.9750-04 

-5.0190-01 

6.5990-01 

-1.3020-02 

-5.1350-04 

-4.9720-02 


1.6390-04 

-4.1080-06 

3.7460-06 

-2.5260-06 

-2.8330-07 

-1.3400-01 

9.9090-01 

-1.0450-02 

-1.0510-03 

-7.1000-03 


1.3490-05 

-3.3180-07 

9.3730-07 

-2.7780-07 

-4.2520-09 

-9.9010-01 

1.3160-01 

-4.4640-03 

7.4310-04 

4.8420-02 


1.7850-05 

-2.4890-08 

4.5080-07 

-2.0760-08 

-3.2890-11 

-9.6490-01 

-5.8210-02 

-4.0560-02 

2.3280-03 

1.0820-02 


3.7710-05 

-1.3620-07 

1.1860-06 

-1.2180-07 

-3.8520-10 

-2.4390-01 

4.9000-02 

-1.5170-02 

-9.7790-04 

-4.1130-02 


-1.7830-05 
6.5150-09 
-3.7750-07 
1. 4390-09 
5.1370-12 
-7. 8680-02 
-5.7340-02 
-2.1290-03 
-3.7840-05 
9. 9520-01 


St2 eigenvalues 


Airframe 


Engine 


Re. Im. 


-3.6740-03 

-1.5510+00 

-1.5510+00 

-1.6720+00 

-1.6720+00 

-9.8330-01 

-3.3740+00 

-1.0070+01 

-1.0070+01 

-2.6240+02 


0.0 

9.1380-02 

-9.1380-02 

1.4500+00 

-1.4500+00 

0.0 

0.0 

1.9850+00 

-1.9850+00 

0.0 


St 2 modal matrix 


1.4820-01 

-7.7140-05 

3.4750-07 

-9.4590-05 

1.2610-03 

5.0720-01 

8.4840-01 

-8.3750-03 

-1.1900-03 

3.0810-02 


-8.5830-03 

-7.8010-04 

7.0120-04 

-4.8510-04 

-5.3170-05 

-1.7900-01 

9.7610-01 

-1.0600-02 

-1.0010-03 

-9.2400-03 


8.5600-03 
8.3670-04 
-9.0810-04 
5.5710-04 
4.4760-05 
3.2740-02 
1 . 1740-01 
-7.7690-04 
-1.5190-04 
1.7650-03 


1.0030-01 

1.5200-02 

-1.6990-02 

1.2400-02 

3.9670-05 

3.5200-01 

4.0130-01 

7.5540-03 

-6.2750-04 

5.4740-02 


8.1780-02 
1. 3770-04 
2.2300-02 
-2.5830-03 
4.6660-04 
3.9610-01 
-7.3220-01 
8.8930-03 
6.5130-04 
3.2640-02 


5.3240-02 

3.5170-03 

-1.4930-03 

1.5180-03 

5.3980-04 

2.4240-01 

9.6860-01 

-1.0610-02 

-1.2460-03 

-3.9250-03 


-2.5210-05 

-2.0890-06 

6.0410-06 

-1.7900-06 

-2.3490-08 

9.9010-01 

-1.3140-01 

4.4550-03 

-7.4420-04 

-4.8460-02 


-2.6290-05 

2.6370-08 

5.0440-08 

1.1980-08 

2.7660-10 

9.9520-01 

4.4460-02 

4.3040-02 

-2.0180-03 

-4.2840-04 


-3.2020-05 

-6.8500-08 

8.9220-07 

-8.6210-08 

5.2140-10 

0.0 

-6. 1780-02 
4.7680-03 
1.5190-03 
4.2550-02 


-1.7830-05 

6.6140-09 

-3.5700-07 

1.3600-09 

5.3170-12 

-7.8680-02 

-5.7340-02 

-2.1290-03 

-3.7840-05 

9.9520-01 


lo 


TABIiE 15. Suboptimal control eigen-structure. 





Airframe 


St4(.l) eigenvalues 

Engine 


Re. 

/ “4.800D-03 
I “1.395D+00 

< -1.395D+00 
I -1.8920+00 

-1.8920+00 
f -7.1050-01 
I -3.5190+00 

< -9.0210+00 
i -9.0210+00 
\ -2.5900+02 


Im. 

0.0 

1 . 1050+00 
-1.1050+00 
1 .0830+00 
-1.0830+00 
0.0 
0.0 

8.0330-01 

-8.0330-01 

0.0 


St4(.l) modal matrix 


2.8230-01 

-1.4440-04 

9.1770-07 

-1.9120-04 

2.5910-03 

4.6600-01 

8.3510-01 

-3,6870-02 

-4.0320-03 

-6.6270-02 


1.5230-01 

8.4930-03 

8.6550-03 

3.5710-03 

7.3120-04 

3.9530-01 

-4.5950-01 

1.7310-03 

-8.6870-04 

2.1800-02 


-7.0520-02 

-1.4020-02 

2.1170-02 

-1.2350-02 

2.5970-04 

1.7660-01 

-7,5460-01 

2.6320-04 

1.0500-03 

-4.8910-02 


1.4570-01 

2.1450-02 

-3.1740-02 

1.7210-02 

2.3610-04 

1.0680-01 

7.3800-02 

6.0840-03 

-9.2410-04 

6.5710-02 


1.0480-01 

2.7500-03 

2.0080-02 

-7.6520-04 

6.2830-04 

4.0280-01 

-8.8440-01 

4.9250-03 

-1.3200-04 

7.2340-03 


8.3380-02 
4.6260-03 
-6. 9660-04 
9.8050-04 
1.3630-03 
6.0210-01 
7.9120-01 
-2.2760-02 
-1.8630-03 
-6.3140-02 


-2.3840-03 

-5.1050-04 

1.5490-03 

-4.4020-04 

-5.2990-06 

9.9140-01 

-1.2560-01 

6.8860-03 

-7.1510-04 

-3.4770-02 


-2.7030-05 

6.0560-06 

-5.2320-05 

6.0720-06 

-7.9190-10 

9.9650-01 

2.6260-02 

3.5410-02 

-1.9960-03 

-2.2250-03 


-4.2400-06 

-2.9420-06 

3.2490-05 

-3.0610-06 

3.4140-09 

-6.5170-02 

-2.1620-02 

-1.9890-03 

5.5900-04 

1.7450-02 


-1.7100-05 

1.6060-08 

1.6250-06 

-6.2760-09 

2.2910-11 

-8.1050-02 

-5.6340-02 

-2.2020-03 

-6.1650-06 

9.9510-01 


St4(.01) eigenvalues 


Airframe 


Engine 


Re. 

-4.7350-03 
-1.4500+00 
-1.4500+00 
-1.8630+00 
-1,8630+00 
-7.6890-01 
-3. 3520+00 
-1.0150+01 
-1.0150+01 
-2.5770+02 


Im. 

0.0 

1 . 0990+00 
-1.0990+00 
1.0520+00 
-1.0520+00 
0.0 
0.0 

1 . 9820+00 
-1.9820+00 
0.0 


St4(.01) modal matrix 


4.2230-01 
-2. 1730-04 
1.3460-06 
-2.8430-04 
3.7580-03 
7.2810-02 
8.9790-01 
-4.1320-02 
-5.2530-03 
-9.1140-02 


1.2380-01 
1.7820-02 
-2.0980-02 
1.4420-02 
-9.6660-06 
-2.4080-02 
6. 1800-01 
-5.5390-03 
-1.4750-03 
2.1060-02 


1.3640-01 
4.6460-03 
1.5770-02 
6.0550-05 
8.3260-04 
3.6970-01 
-6.6720-01 
1 .1880-03 
-4.6230-04 
-4.3570-03 


1.4630-01 

2.1690-02 

-3.2580-02 

1.7450-02 

2.2470-04 

1.2760-01 

2,0480-01 

-2.1390-03 

-1.2570-03 

2.0560-02 


1.1480-01 

3.9060-03 

1.8230-02 

6.9920-05 

6.7360-04 

3.5910-01 

-8.8070-01 

4.7320-03 

-3.8430-05 

4.0670-05 


1.3440-01 
7.7330-03 
-1.6170-03 
2.1030-03 
1.9450-03 
4.4000-01 
8.8580-01 
-2. 2500-02 
-2.2250-03 
-5.5690-02 


3.8840-03 

7.8610-04 

-2.2480-03 

6.7060-04 

9.1520-06 

-9.8900-01 

1.3950-01 

-4.6490-03 

6,6190-04 

4.8700-02 


-2.0740-05 
-3.5290-06 
4.4530-05 
-3.7510-06 
5.5550-09 
9. 9610-01 
3.4770-02 
4.2190-02 
-1.9930-03 
-3.8410-03 


-2.9690-05 
-3.2560-06 
2.5670-05 
-3.2620-06 
1. 2400-09 
r3. 0120-04 
-5.6470-02 
3.9340-03 
1.4250-03 
4.0300-02 


-1.7100-05 

1.6120-08 

1.6430-06 

-6.3780-09 

2.3180-11 

-8.1400-02 

-5.6370-02 

-2.2270-03 

-6.2120-061 

9.95ID-OI! 


TABLE 15. concluded 



Airframe 

St3(.l) eigenvalues 

Engine 


Re. 

! -4.7*560-03 
-1. 3220+00 
-1 .3220+00 
-2.0020+00 
-2.0020+00 
/ -7.5920-01 
\ -3.3310+00 
< -1.0080+01 
i -1.0080+01 
^ -2.6220+02 


Im. 

0.0 

1.0520+00 

-1.0520+00 

1.1970+00 

-1.1970+00 

0.0 

0.0 

1.9760+00 

-1.9760+00 

0.0 


St3(.l) modal matrix 


4. 4710-01 
-2.2970-04 
1.4340-06 
-3.0150-04 

4.0080-03 
3.3540-02 
8,8750-01 
-4.4370-02 
-5.6890-03 
-9.6220-02 


1.4070-01 
6.8570-03 
9.1090-03 
2.2200-03 
7.6330-04 
1 .7210-01 
-6.4640-01 
1.4580-03 
-3.8150-04 
8. 0760-03 


-7.9450-02 
-1,3430-02 
1.7480-02 
-1.1450-02 
2.1110-04 
-9.2820-02 
-7. 1850-01 
8.2060-04 

7.0080-04 
-2.7770-02 


1.2150-01 
1.9910-02 
-3.2900-02 
1.6420-02 
1.4830-04 
3.4550-01 
-1 .2060-01 
9.6770-03 
1.7470-04 
4.0970-02 


1.0080-01 

3.3250-03 

1,9620-02 

2.0800-05 

5.2700-04 

-2.4350-02 

-9.1460-01 

3.1260-03 

1.7660-04 

1.1740-02 


1.4180-01 
8.1150-03 
-1.6170-03 
2.1310-03 
2.0930-03 
4.0210-01 
9.0270-01 
-2.2560-02 
-2.2810-03 
-5. 2800-02 


-1.7080-03 

-3.4220-04 

9.7130-04 

-2.9150-04 

-4.0410-06 

9.9060-01 

-1.2700-01 

3.8150-03 

-7.5590-04 

-5.0330-02 


-4.6990-05 

1.5980-05 

-1.5860-04 

1.6360-05 

-1.2460-08 

9.8650-01 

3.2160-02 

4.4210-02 

-1.7240-03 

8.2450-03 


-6.6430-06 

-2.6320-06 

6.4680-05 

-3.2100-06 

1.2780-08 

-1.3400-01 

-6.5270-02 

-1.6850-03 

1.7200-03 

3.9950-02 


-1.7070-05 

1.6970-08 

1.7940-06 

-6.8430-09 

2.4110-11 

-7.9070-02 

-5.7330-02 

-2.2400-03 

-3.7830-05 

9.9520-01 


Airframe 

St3(.01) eigenvalues 

Engine 


-4.7270-03 
-1.4470+00 
-1.4470+00 
-1 .8810+00 
-1.8810+00 
-7.7130-01 
-3.3770+00 
-1.0070+01 
-1.0070+01 
-2. 6240+02 


0.0 

1.0800+00 

-1.0800+00 

1.0890+00 

-1.0890+00 

0.0 

0.0 

1.9840+00 

-1.9840+00 

0.0 


St3(.01) modal matrix 


4.4940-01 

-2.3130-04 

1.4290-06 

-3.0240-04 

3.9920-03 

-1.9930-02 

8.8770-01 

-4.2430-02 

-5.6370-03 

-8.8400-02 


1.6820-01 
9.6530-03 
8.9240-03 
4.2060-03 
8. 2020-04 
2.7600-01 
-5.2190-01 
-7.3610-04 
-8.4860-04 
3.5780-03 


-8.2410-02 

-1.5920-02 

2.4660-02 

-1.3900-02 

2.4080-04 

6.9600-02 

-7.8090-01 

4.3220-03 

1.1870-03 

-2.2400-02 


1.5270-01 
2. 2680-02 
-3.3720-02 
1.8250-02 
2.3380-04 
1.5630-01 
1.1700-01 
1.3780-04 
-1.0290-03 
2.4360-02 


1.1470-01 

3.2890-03 

2.0900-02 

-5.4660-04 

6.7690-04 

2.4280-01 

-9.2920-01 

3.4740-03 

-5.6340-05 

1.0990-03 


1.4720-01 

8.4870-03 

-1.7950-03 

2.3270-03 

2.1210-03 

3.8760-01 

9.0800-01 

-2.3460-02 

-2.4010-03 

-5.4400-02 


-2.2660-03 

-4.6240-04 

1.3340-03 

-3.9510-04 

-5.2900-06 

9.9110-01 

-1.2340-01 

4.3110-03 

-7.4850-04 

-4.9270-02 


-2.7290-05 

-7.9520-07 

1.3650-05 

-9.0250-07 

2.4480-09 

9.8550-01 

3.5360-02 

4.3280-02 

-1.7860-03 

5.5400-03 


-2.6630-05 

-2.2430-06 

2.1360-05 

-2.2980-06 

1.9280-09 

-1.3880-01 

-6.7360-02 

-1.2840-03 

1.7850-03 

4.2170-02 


-1,7680-05 

8.8150-09 

9.8390-08 

-3.7490-10 

9.3000-12 

-7.8680-02 

-5.7340-02 

-2.1290-03 

-3.7840-05 

9.9520-01 


in 


TABLE 15. continued 
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FIGURE 2. Separate and integrated system transient response. 
Offset in angle of attack of 0.0008 rad. 
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FIGURE 3 



Separate and integrated system transient response 
Offset in altitude of 0.137 km. 
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FIGURE 15, concluded 
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